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ABSTRACT 


This research has examined the benefits of using the Robust Symmetrical Number 
System (RSNS) to resolve ambiguities in phase sampling interferometry. A compact, 
high resolution direction finding antenna architecture based on the RSNS was developed 
to demonstrate experimentally the elimination of phase errors using a minimum amount 
of hardware. Previous work has determined that phase errors in the system will degrade 
the system performance. Several improvements were made to the original RSNS 
prototype antenna to provide enhanced performance. Adding isolators and supplementing 
the ground plane with copper tape (between the antenna elements), a reduction in the 
mutual coupling effects was accomplished. Mounting the microwave components on a 
brass plate also reduced errors contributed by vibrations and temperature. Tailor cutting 
all semi-rigid coaxial lines also helped reduce the number of connectors required to 
assemble the microwave circuit, also a source of phase errors. Matching the front-end 
amplifiers in each amplification stage rather than matching the characteristics of two 
cascaded amplifiers in each signal line has reduced relative phase errors between 
channels as well as matching the power outputs of the amplifiers. Two printed circuit 
boards were designed and built for the RSNS signal processor. The printed circuit boards 
provide a decrease in the electrical noise floor over the original design (assembled on 
breadboards). The new design has reduced the phase errors that were present in the first 
prototype system. The RSNS signal processing technique is able to provide a high- 
resolution phase sampled direction finding capability with an angular resolution of 1.9 


degrees by using only three receiving elements (two interferometers). 
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I. INTRODUCTION 


A. DIRECTION FINDING ANTENNAS 


Direction Finding (DF) antennas are used by many people. Law enforcement, 
military, wildlife managers and telecommunications personnel use DF systems for 
locating and tracking different types of emitters. A wide variety of techniques are used to 
accomplish the DF task. Spinning DF, time difference of arrival (TDOA), and phase 
sampled interferometry are just a few examples that are commonly used. All these 
systems are passive; meaning they do not generate their own signal to locate a target. 
Therefore most of these systems will only return a line of bearing to a transmission. 

Phase sampled interferometry is a very attractive technique since a direction of 
arrival can be obtained by the phase information stored in a plane wave. Amplitude 
information is not required. Within the phase sampled interferometry field of direction 
finding, several new techniques have been investigated. 

The ESPIRIT-based two-dimensional arrival estimation scheme achieves aperture 
extension (interferometry baseline extension) using a sparse uniform rectangular array of 
electromagnetic vector sensors spaced much farther apart than a half-wavelength [1]. An 
electromagnetic vector sensor is composed of six spatially co-located, orthogonally 
orientated, diversely polarized antennas, distinctly measuring all six electromagnetic-field 
components of an incident multisource wavefield. The direction of arrival from each 
incident source is estimated from the source's electromagnetic-field vector component 
and serves as a course reference to remove the cyclic phase ambiguities in ESPIRIT's 
eigenvalues when the intervector sensor spacing exceeds a half-wavelength. 

A polynomial rooting approach to super-resolution array design is concerned 
with the design of an array that satisfies prespecified performance levels, such as 
detection-resolution thresholds and Cramér-Rao bounds on error variance [2]. The sensor 
placement problem is formulated in the framework of subspace-based DF techniques and 


a novel polynomial rooting approach to the design problem, based on the new concept of 
] 


the "sensor locator polynomial." This polynomial is constructed using prespecified 
performance levels, and its roots yield the sensor locations of the desired array. The 
distinguishing feature of this technique is that it hinges on the properties of the array 
manifold. 

The extended phase interferometry technique [3] incorporates both calibrated 
phase and amplitude response data from the antenna arrays. This technique appropriately 
weights the square of the baseline phase differences with the antenna gains. The 
incorporation of amplitude data provides significant performance improvement over 
phase-only interferometry; however, it requires a modest increase in computational 
complexity. 

This thesis research improves on a phase sampling interferometer approach that 
can be easily incorporated into the established techniques to provide a high-resolution, 
small baseline array with fewer number of phase sampling comparators [4]. The approach 
is based on preprocessing the received signal using the robust symmetncal number 
system (RSNS). The RSNS preprocessing is used to decompose the spatial filtering 
operation into a number of parallel sub-operations (moduli) that are of smaller 
computational complexity. Each sub-operation 1s a separately configured interferometer 
that symmetrically folds the detected phase difference with folding period equal to 2Nm; 
where N is the number of interferometers that are used within the linear array. A small 
comparator ladder mid-level quantizes each folded phase response. Consequently, each 
sub-operation only requires a precision in accordance with that modulus. A much higher 
DF spatial resolution is achieved after the N different moduli are used and the results of 
these low precision sub-operations are recombined. By incorporating RSNS 
preprocessing concept, the field of view of a specific configuration of interferometers and 
phase sampling comparator arrangements can be analyzed exactly. Experimental results 


for an improved 6-bit RSNS array are presented. 


iB: PRINCIPAL CONTRIBUTIONS 


The focus of this thesis research has been to improve the design and performance 
of a previous RSNS direction finding array. Microwave component performance was re- 
evaluated and new features were introduced. A new RSNS DF system was fabricated 
and tested. 

Initial design efforts focused on the microwave receiver. The previous 
prototype's performance was marginal. Its design injected distortions into the folding 
waveforms, causing angle of arrival reporting errors. The existing component 
performances were re-examined and a new design was constructed. 

The previous RSNS digital processing circuit was designed and constructed 
without its performance ever being validated. Computer simulations of the new circuit 
were conducted to evaluate the current design. During the simulations, it was discovered 
that by adding a timing circuit to clock in comparator outputs to the RSNS encoders, the 
performance of the system could be improved. 

The new RSNS processing circuit was laid out for construction on a printed 
circuit board. The circuit board was milled on campus in the Physics department 
electronics workshop. 

The improved prototype was tested. Data collection was conducted by collecting 
mixer output voltages and RSNS processor outputs. The measured patterns were used to 
study the effects of the new design features on mutual coupling, amplifier harmonics, 


temperature, vibrations and power levels. 


Cc. THESIS OUTLINE 


This thesis presents the theoretical background and design equations needed to 
construct a phase sampled interferometer, as well as previous test arrays that have been 
built and tested. A new array was designed to eliminate phase errors that were present in 


the preceding systems. 


Chapter II is a review of phase sampled interferometry and the prototype antenna 
design elements. Chapter Il introduces the Optimum Symmetrical Number System 
(OSNS), taken from the work of Thomas Hatziathanasiou [4]. It presents the design of a 
DF system based on the OSNS. Chapter IV, taken from the work of David Wickersham 
[S], gives an overview of previous work done with the Robust Symmetrical Numbering 
Systems and its application in DF antenna systems. Chapter V presents my design data 
for the antenna and RF receiver components. Chapter VI is a detailed description of my 
RSNS signal processing design. Chapter VII is an explanation of the testing procedures 
for the antenna and a presentation of the experimental results. Chapter VIII gives some 


concluding remarks and summarizes the next step that is needed. 


Il. PHASE SAMPLED INTEFEROMETRY 


A two-element linear interferometer is shown in Figure 2.1. The two antenna 
elements are spaced a distance d apart and the incident plane wave arrives with bearing 


angle 8,[5]. In this phase monopulse configuration the angle 9, is measured from the 
perpendicular to the baseline axis and can take on values 1/2 >6, >-1/2. The phase 


difference between the two elements is 
ph) 
Ay=VW,-VW, Sea OU, (2.1) 


where d is the element spacing and A is the wavelength, and is a function of the incidence 


angle of the wave. 
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Figure 2.1: Two Element Interferometry. 
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The signals received by the antenna elements are mixed (multiplied together and 
lowpass filtered) resulting in an output signal whose frequency is the difference of the 
two input signal frequencies. Because the signals have the same frequency, the mixer 
output voltage is a value proportional to the difference between signal phases wy, and 
y,. This difference is not purely Aw because of the time delays due to the different 


transmission line lengths, oj, from each antenna element. However, these line lengths are 
known and can be compensated for in the angle estimate. Let the signals from the two 


antenna elements be 


v,(t) = Vcos[2nft + 0, (t)] (2.2) 


and 
vV,(t) = Vcos[2nft + 6, (t)], (2.3) 


where V is the maximum value of the voltage at the antenna elements, t is time and f is 


frequency. Let w, be the phase difference between the transmission lines to the two 


elements. The lowpass mixer output voltage is 
WE i irda 
Vor = "ND = ace ag ey teal (2.4) 


which contains the plane wave angle of arrival (AOA) information. For values of 
d=A/2 and y, =0, Ad=m7sin(@,). As the AOA @, varies from —1/2 to 1/2, the 
phase difference Advaries from —7 to % as shown in Figure 2.2. The output voltage 
from the phase detector is also a function of the phase difference and is a symmetrical 
folding periodic waveform. Together these relationships give the phase detector output 


voltage as a symmetrical function of the AOA as shown in Figure 2.3 for d=A/2 [6]. 
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Figure 2.2: Mixer Output Voltage Versus Phase Difference, d=A /2 [from 6]. 
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Figure 2.3: Mixer Output Voltage Versus Angle of Incidence, d=/ /2 [from 6]. 


Ambiguities are generated for baselines where d>A/2. That is, the phase 
detector output voltage is highly ambiguous with a single phase corresponding to many 
angles of arrival. The number of folding periods n that occur within an AOA of 7 radians 


1S 


a (275) 


For example, with d = 7.54, n=15 folds are available as shown in Figure 2.4. The folding 
period is not constant but grows larger in proportion to the angle off of broadside because 


of the sin(@, ) dependence in (2.4). 
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Figure 2.4: Mixer Output Voltage Versus Angle of Arrival for d = 7.5A [from 6]. 


The ambiguities within the symmetrical folding waveforms represent the phase difference 
between the elements and can be resolved by using additional interferometers in the 
linear array. Typically, each interferometer in the linear array symmetrically folds the 
phase response with the folding period between interferometers being a successive factor 


of two, or d, =2d, =4d, =8d,, where d, represents the element spacing. High-speed 


binary comparators are used to produce a digital output. The folding waveforms are 
shifted appropriately using a phase shifter in each channel to achieve a Gray code result. 
The folded output from each phase detector is then quantized with a single comparator 
with a normalized threshold level T=0.5. Together, the comparator outputs directly 
encode the signal's AOA in the Gray code format. This approach makes use of the 
periodic dependence of the interferometer's phase response on the applied plane wave's 


AOA and the distance between the elements of each interferometer. One of the major 
8 


limitations associated with this approach is the achievable resolution. For the folding 
periods to be a successive factor of 2, the distance between the elements must also be 
doubled. That is, an 8-bit DF antenna using the previous scheme would require element 
spacings 4/4,A/2,A,2A,42...,32A with a total baseline length of 32A. This distance- 
doubling of the element spacings requires complex analog hardware, adversely affects the 
physical implementation of the DF architecture and ultimately constrains the achievable 


resolution. 
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IW. OPTIMUM SYMMETRICAL NUMBER SYSTEM DIRECTION 
FINDING ANTENNA 


A. OPTIMUM SYMMETRICAL NUMBER SYSTEM 


Chapter IJ summanzes Thomas Hatziathanasiou's work in symmetrical number 
systems [4]. The optimum symmetrical number system (OSNS) is composed of a 
number of pairwise relatively prime (PRP) moduli m,. The integers within each OSNS 
modulus are representative of a symmetnically folded waveform with the period of the 
waveform equal to twice the modulus, t.e., 2m,. For m given, the integer values within 
twice the modulus are given by the row vector [6, 7] 


x. =[0J,---,;m—1,m-1,---,1,0] for m=m;, i=1, 2, ... (3.1) 


m 


Figure 3.1 shows part of the OSNS folding waveforms and the integer values within the 
modulus for m,;=5 and m2=6. The horizontal axis represents the normalized input. The 
vertical lines represent folding waveform reference levels. The numbers at the top of the 
figure represent the number of reference levels that are crossed by the folding waveform 


for a given input value. The period of one complete fold is equal to 2m, =10 and 


2m, =12. 
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Figure 3.1: Folding Waveforms and Integer Values Within Moduli m, = 5 and m, = 6 for 
the OSNS [from 5]. 


Due to the presence of ambiguities, the integers within the vector derived by equation 
(3.1) do not form a complete system of length 2m by themselves. The ambiguities that 
arise within the modulus are resolved by considering the paired values from all channels 
together. If N pairwise relatively prime moduli are used, the dynamic range (the number 


of unambiguous vectors) of this scheme is 


M=][m,. | (3.2) 


i=] 


This dynamic range is also the position of the first repetitive moduli vector. As shown in 
Figure 3.1, the dynamic range for the m, =5 and m, =6case is M=30. That is, no set 


appears twice for a normalized input in the range (0,29) [6, 7]. 


B. OSNS PROTOTYPE DF ANTENNA 


Figure 3.2 shows the schematic diagram of a N=2 channel OSNS DF antenna. 


The array consists of three elements, with all channels sharing a common element. With 


12 


the proper distances between elements, each interferometer folds the phase response at 
2m;. The phase response from each channel is quantized using mj-1 comparators. The 
integers shown in Figure 3.1 represent the number of comparators ON due to the phase 


voltage exceeding the comparator matching threshold voltage. 
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Figure 3.2: Optimum Symmetrical Number System Antenna Architecture [from 4, 8]. 


The distance between each pair of elements is a function of the modulus, m;, The phase 
response must have the correct number of folds to cover the dynamic range. For a given 


modulus m;, the number of folds within the dynamic range is 


nei / 2m (3.3) 
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Since each fold of the array corresponds to a change in phase of m for d=A/2, the 
required spacing between the reference element and the respective element of channel i is 


given by [5] 


d.=n.—= 


1 i 


(3.4) 


wm | > 


MA 
4m, 


The bin widths are uniform in size in sin(®,) space with a width 2/M. Therefore, the 


spatial resolution is given by 





oer 9 Ik - 
—— arcs EMS — resi & ag, (3.4) 


where ke {0,1,---M—1} and k=0 is the bin that starts at 6, =—-1m/2. The variable€is a 


scale factor that can be used to compress the encoded field of view symmetrically about 
broadside. The phase waveforms are not ideal over 180 degrees because of the element 
pattern degradation at wide angles. Since the waveforms at wide angles are not usable, 
the SNS can be remapped to the usable portion of the waveform to increase resolution. 


The scale factor is defined to be 


_d__ sin(Q) 


d’ sin(@’)’ 





(3.5) 


where @ corresponds to phase angle for the unscaled configuration, and @ is the desired 
phase angle for the compressed configuration. Note that the spatial resolution is not 
uniform over the range of -1/2 to 7/2. The resolution is finer than r at broadside (6 = 0°) 


and increases with angle off broadside, because of the sin(®, ) dependence in (2.4). 


Cc. EXPERIMENTAL RESULTS 
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A linear array based on the moduli m; = 6 and m2 = 11 was designed, fabricated 
and tested at a frequency 8.5 GHz [4, 6]. The radiating elements are open-ended 
waveguides as shown in Figure 3.2. For moduli m; = 6 and m2 = I1, d)=3.82 in and 
d,=2.08 in. With M=66, the number of folding periods is n;=5.5 and n2=3. To provide an 
adequate signal-to-noise ratio, a low-noise amplifier is included at the output of each 
interferometer element. Since the common element splits the signal into N paths, an 
attenuator is placed in the other branches to balance the amplitudes. The amplifiers 
operate in saturation so that the mixer input signal levels are independent of the angle of 
incidence. A fixed phase shifter is also included in one branch of each interferometer so 
that the symmetrically folded phase response waveforms from each mixer may be 
aligned. This alignment insures that the comparators in the digital processor properly 


sample the phase waveform and encode it in the OSNS. 


The anechoic chamber facility at the Naval Postgraduate School was used for the 
antenna pattern measurements for each pair of elements. The measured mixer outputs are 
shown in Figure 3.3 [4, 8]. Since the array is not useful in the endfire regions, the phase 
response from each channel was minimized at an input AOA of -50°. This is the point 
where the digital code starts. The locations of the voltage minima are controlled by 
introducing the appropriate phase shifts to the signals at a point between the elements and 
the mixer. A simulation of the transfer function for the OSNS array is shown in Figure 


3.4 for €=1. The measured transfer function (output of digital hardware) is shown in 


Figure 3.5 (4, 8]. 
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Figure 3. 3: Measured Mixer Output Voltages for OSNS Array [from 4, 8]. 
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Figure 3. 4: Simulated Transfer Function Using the OSNS Measured Mixer Output 
[from 4, 8]. 
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Figure 3. 5: Measured Transfer Function of OSNS Prototype Antenna [from 4, 8]. 


The spikes in the transfer function are AOA encoding errors. The phase waveforms must 
cross the comparator thresholds (code transitions) at the same time for all channels. The 
OSNS can present a problem at each of these specific AOAs. When some comparators, at 
positions of change, do change, while others do not, the recovered amplitude will have a 
large error. Encoding errors can also occur due to imperfections in the transmission line 
phases. These errors however, can be isolated with additional comparators and corrected 
using interpolation circuits as described in [6]. An alternate approach, however, is to 
develop a symmetrical number system that eliminates the possibility of this type of 
encoding error. That is, a symmetrical number system is needed where only one 


comparator threshold is crossed at any one particular code transition. 
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IV. ROBUST SYMMETRICAL NUMBER SYSTEM 
DIRECTION FINDING ANTENNA 


A. ROBUST SYMMETRICAL NUMBERING SYSTEM 


Chapter [IV summarizes David Wickersham's work in symmetrical number 


systems [5]. In the RSNS, N different periodic symmetrical waveforms are used with 


pairwise relatively prime integers m;, m2, ...,. My. The RSNS is based on the following 
sequence 
X_, = (0, 1, 2,---, m-1, m, m-1,---, 2, 1], (4.1) 


where X_, 1S a row vector and m is a positive integer (m > QO) [9]. In an N-channel 


RSNS, where N = 2, the basic sequence for the i” channel (modulus m) is 
Xe = [o, O,---, CQ. 0, iP l, ney l, l, mare 00 m,°°*,™M, mM, res Ie Poo. l, 1]. (4.2) 


In this sequence each value in the x,, row vector is put N times in succession. This 


sequence 1s repeated in both directions, forming a periodic sequence with the period 
Prsns = 2mN. (4.3) 


Considering a single channel, the discrete states of the robust symmetrical number system 


can be expressed as [9] 





a Son = Nims. +1 
g = (4.4) 
Am Sele ene ee 
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where g is the n™ term of channel i, m; is the channel modulus, s; is a corresponding 
sequence shift s; = 0, 1, 2,..., N-1 (mod N) and N 2 2 is the number of channels in the 
system. The values {s;, S2,..., Sn} must form a complete residue system modulo N. 
Because of the relative property of the shifts, one of the shift values will be set equal to 0. 
The index n corresponds to the input value. The discrete states of the RSNS are indexed 
using the s; =Q row vector with the index starting from the first zero. Note that the 
largest integer within each periodic sequence is the modulus m. 

An N-channel RSNS is formed of vectors by picking N moduli mj, and N shift 
values sj, 1 <1 S N. Since the fundamental period for channel 1 is 2Nmj, it follows that 
the period for the RSNS vectors must be a multiple of 2Nm;. Therefore the fundamental 
period for the RSNS is [9] 


PFersns = [2m,N, 2moN,..., 2MnN], (4.5) 
where [a}, a2, ..., an] is the least common multiple of aj, a2, .... an. From number theory 
PFrsns = 2N{m1, mo,..., mn]. (4.6) 


To illustrate the RSNS, Figure 4.1 shows the folding waveform and integer values within 
each modulus for m;=5 (shift s;=0) m2=6 (shift s2=1). Note the Gray-code properties 
from one code transition to the next. The thresholds shown on the vertical axis represent 
the integer values within each RSNS modulus. The integer values occur N=2 times in 


succession. 
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Figure 4.1: Folding Waveforms and Integer Values Within Moduli m, = 5 and m, = 6 for 


the RSNS. 


The system dynamic range (M,.,.) of the RSNS is the maximum number of 


RSNS 


distinct vectors without a redundancy. The selection of the shift and the permutations 
among the 2 moduli has no effect on Mis, However, the point indices corresponding 
to M ass (beginning point and ending point) are different. When the channel moduli are 
of the form m, = 2“ -1,m, = 2* im, = 2" +1, the closed-form expression for Masys is 


conjectured to be 


“ 3 15 
Masns = 7m coh +7, (4.7) 


where m; 2 3 [9]. For two channels, with the moduli separated by 3 or more, 


Mass =4m, +2m, -2. (4.8) 
Checking against computer results when m)=8 and m2=17, M = 64, which fits (4.8). It 
can be verified that Meee =27<M =30 (3.2) for the two channel example shown in 


Figure 4.1. 


| 


B. RSNS PROTOTYPE ANTENNA 


To demonstrate the efficiency of the RSNS preprocessing, the design of a 6-bit 
antenna using m, = 8 and m2 = 17 is considered [5]. A schematic diagram of the RSNS 


antenna is shown in Figure 4.2 for shifts of s;=0, s2=1 resolution bins. 
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Figure 4.2: Block Diagram of the Unscaled RSNS Antenna [from 5]. 


The number of folds within each modulus is 


(4.9) 





or nj=2 and n2=0.94. The required distance between antenna elements is 


G09 


ene) (4.10) 
2 4m,NE 





With a scale factor of €= V¥3/2, d\=1.70 in and d)=0.80 in. For a normalized folding 


waveform amplitude (-1 to 1), the threshold for the jth comparator for the modulus m, 


channel is [5] 


V... =cos) ———*z |. (4.11) 


The input signal is applied in parallel to both interferometers. The output from each phase 
detecting mixer is amplitude analyzed using m,; comparators. For this design, there are a 
total of 25 comparators with a maximum of 17 loaded in parallel. An RSNS-to-binary 
logic block translates the RSNS residues (number of comparators ON in each 
thermometer code) to a more convenient representation. The logic block and 


corresponding bin number are shown in Figure 4.3. 
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Figure 4.3: RSNS Logic Block [from 5]. 


A consequence of quantizing the angle of arrival into a bin is that the RSNS 


aA 


system reports any signal that falls within a bin as if it arrived at the bin center, 6, . That 


is to say, for one angle within each bin, the estimation is exact, but for the remaining 
angles a reporting error exists. The AOA resolution or bin width in @-space for the kth 


bin is 


I, = arcsin| €—————_ |—-arcsin| € (4.12) 
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Another advantage of the RSNS is that small phase errors can be tolerated without 
a serious error to a reported AOA. These phase errors cause the folded waveform from 
the mixer to shift, thus changing the timing of the comparator state changes. Because the 
waveforms from the mixers are quantized, the bin width or r, provides some free-play in 
the antenna circuitry and limits AOA reporting errors to the adjacent quantization levels 
(bins). However, there is a limit to this feature. If the phase error becomes too great, the 
comparators will not change state in the correct sequence, which may or may not 
correspond to a code that is mapped in the dynamic range. These errors are similar to the 
errors that the OSNS experiences when the comparators of both channels do not change 
state at the same time. For the system designed, a theoretical phase error of 6° can be 
tolerated in the modulus 17 channel and a phase error of 12° in the modulus 8 channel. 
Figures 4.4, 4.5 and 4.6 show the simulated transfer function. Figure 4.4 is the transfer 
function with no phase error in both channels. Figures 4.5 and 4.6 are the modulus 8 
channel with zero degrees of phase error and 3° and 6° phase error in the modulus 17 
channel. Note how the transfer function becomes jagged and erratic as the phase error is 


increased until it becomes discontinuous at 6°. 
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Figure 4.4: Simulated Transfer Function with No Phase Errors in Either Channel 


[from 5]. 
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Angle of Arrival (degrees) 


17 Channel with 3° Phase Error Introduced [from 5]. 


Figure 4.5: Modulus 
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Figure 4.6: Modulus 17 Channel with 6° Phase Error Introduced [from 5]. 


c. EXPERIMENTAL RESULTS 


A RSNS array based on the moduli m,; = 8 and m2 = 17 was designed, fabricated 
and tested at a frequency of 8.0 GHz [5]. The schematic diagram of the two-channel 
array is shown in Figure 4.2. The radiating elements are printed circuit dipoles. For 


moduli m,; = 8 and m2 = 17, the unscaled distances, € = 1, are d)=1.48 in. and d,=0.70 in. 


With M = 64, the number of folding periods nj=2 and n2=0.94. To provide an adequate 
signal-to-noise ratio, a low-noise amplifier is included at the output of each 
interferometer element. Also included in the RSNS array is a lowpass filter in each 
channel to eliminate the amplifier harmonics that are present. Since the common element 
splits the signal into N paths, an attenuator is placed in the other branches to balance the 
amplitudes. A fixed phase shifter is also included in one branch of each interferometer so 
that the symmetrically folded phase response waveforms from each mixer may be 


aligned. This alignment insures that the comparators in the digital processor properly 


Zi 


sample the phase waveform and encode it in the RSNS. For the m, = 8 channel, 8 
comparators are required. For the m2 = 17 channel 17 comparators are required. The 
EEPROM that was used accommodates 15 inputs. The modulus 17 channel processing is 
conducted in two stages. In the first stage, the 4-bit thermometer to binary encoder maps 
the 15 least significant codes out of the 17 comparators (ones with the smallest 
threshold). The remaining two comparators are mapped on the second EEPROM along 


with the 4 bit binary encoder output and all modulus 8 comparators. The digital processor 


puts out a bin number k corresponding to the estimated AOA 6. The estimated AOA for 


bin k 1s given as 
6, = sin (E21 -€). (4.13) 


Broadside incidence for the RSNS antenna corresponds to the transition between Bins 31 
and 32. For the modulus 17 channel at slightly negative angles (bin 31), the ninth 
comparator is ON and for slightly positive angles (bin 32), the ninth comparator is OFF. 
That is, the phase value at 8, =0 corresponds exactly to the threshold of the ninth 
comparator. From (4.11) the phase corresponding to the threshold value Vj79 is 1/2. 


This is the phase value required for the modulus 17 channel at 8, =Qand is obtained 
using the phase adjuster shown in Figure 8. In the modulus 8 channel, 8, = Ooccurs 


midway between comparator 7 and 6 thresholds. Since V,,=cos”’ es) and 


Ve Sicos (22), the required phase value is <2 = 4 and is also obtained using the modulus 


8 phase adjuster. 
The normalized mixer folding waveform outputs are shown with the simulated 


waveforms in Figures 4.7 for the m)=8 channel and 4.8 for m2=17 channel for§ =1.0. 


The basic features correspond to predicted curves. The measured folding outputs from 
the shift and bias amplifiers in both channels are shown in Figure 4.9. The digital circuit 
maps the phase difference at the mixer output into an estimated angle of arrival. The 


predicted output of the unscaled prototype array is reproduced as Figure 4.10. The 
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measured mixer output contains phase errors, which result in differences between the 


simulated and predicted data. These phase errors cause the thermometer code to 


incorrectly map the incident wave to the proper angle of arrival. The measured transfer 


function of the unscaled prototype antenna is shown in Figure 4.11. 


Figure 4.7: 


Normalized Mixer Output 
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Angle of Arrival 


Normalized Mixer Output and Simulation Waveforms for RSNS Array with 


m,=8 [from 5]. 
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Figure 4.8: Normalized Mixer Output and Simulation Waveforms for RSNS Array with 
m,=17 [from 5]. 
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Figure 4.9: Measured Folding Waveform Outputs from the Shift and Bias Amplifiers 
[from 5]. 
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Figure 4.10: Predicted RSNS Transfer Function [from 5]. 
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Figure 4.11: Measured RSNS Transfer Function [from 5]. 


The large discontinuities in the transfer function are due to phase and amplitude 


errors, which arise from several sources. Mutual coupling 1s a common cause of phase 
3] 


front distortion in small arrays of closely spaced elements. In order to increase the 


distance between elements and thus reduce the mutual coupling, a scale factor € can be 


introduced. The scale factor is the ratio of the element spacings of the scaled and 
unscaled arrays. Reducing the scale factor narrows the mapable field of view of the 


antenna and increases the resolution within the mapable field of view. For a scale factor 
of E= 3 /2, the maximum mapable aperture (6,,,) 1s 60°. The predicted transfer 


function for the scaled array is shown in Figure 4.12. The corresponding measured 
transfer function is reproduced as Figure 4.13. The mixer outputs contain significant 
phase errors, and the increased resolution of the antenna limits the ability to process the 
errors. The steeper slope of the scaled transfer function is clearly visible when 


comparing 4.10 and 4.12. 


Reported Angle of Arrival (degrees) 





-80 -60 -40 -20 0 20 40 60 80 
Angle of Arrival (degrees) 


Figure 4.12: Predicted RSNS Transfer Function for Scale Factor € = V3/2 [from 5]. 
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Figure 4.13: Measured RSNS Transfer Function for Scale Factor € = V3/2 [from 5]. 


The transfer function for the scaled array of dipoles does not show any significant 
improvement due to the reduction of mutual coupling. To verify this, a second array 
comprised of open-ended waveguide elements was constructed and tested. These 
waveguide elements have significantly less H-plane mutual coupling than the dipoles, but 
the transfer function showed no improvement over with printed circuit (dipole) elements. 
This indicates that mutual coupling is not the primary cause of the phase errors. The 
microwave circuit has three likely sources of errors: (1) the low noise amplifiers, (2) 
coaxial cables and phase adjusters, and (3) the rigid connectors. The prototype uses two 
cascaded LNAs operating in saturation in each channel. Harmonics generated by the first 
amplifier are amplified by the second amplifier. The second amplifier also introduces its 
own harmonics. Filters were used to suppress harmonics. 

The coaxial cables and phase adjusters are another source of error within microwave 
circuits. The ngid coaxial cables are coarsely trimmed at 8.0 GHz, and the phase shifters 
are used for fine trimming. The cable lengths are long and subject to thermal expansion 


and contraction, particularly in the vicinity of the amplifiers. Finally, all of the devices in 
a3 


the circuit are jorned with SMA type connectors. The circuit was repeatedly assembled 
and disassembled for troubleshooting, and consequently the quality of the connections 
were degraded. Phase errors of + 2° can easily be introduced at each connection. Thus it 
appears that many small errors accumulate in a manner such that the 6° and 12° error 


thresholds are exceeded, which results in a degraded transfer function. 
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V. IMPROVED RSNS RF SYSTEM DESIGN 


Previous research has introduced the Optimal Symmetrical Number System and 
the Robust Symmetrical Number System and shown that they are able to efficiently 
resolve ambiguities in phase sampling interferometry. The prototype arrays that were 
built, however, both showed that significant errors are introduced due to the nature of the 
components and construction of the system. Non-linear components, poor connections 
and mutual coupling are prime contributors to the phase errors present in the folding 
waveforms. This thesis research takes a step forward to validate the RSNS design and to 
eliminate the phase errors in the system. 

The RF circuit shown in Figure 5.1] receives a transmitted signal with a set of 
micro-strip antennas (see Figure 4.2). The signal is passed through a set of isolators and 
then to the first amplification stage. Band-pass filters are used to remove harmonics and 
the signal is then sent to the second amplification stage. Low-pass filters are used to 
remove harmonics from the second stage amplifiers. A power splitter is used to divide 
the reference signal that is passed on to the RF side of the balanced mixers. The signals 
from each channel are passed through phase shifters and then to the local oscillator (LO) 
side of the mixers. Since the frequency at the RF and LO inputs are the same, the mixer 
output is a DC value that is a function of the AOA. This DC function is referred to as the 
folding waveform that is amplitude analyzed by the RSNS signal processor in order to 


determine the AOA. 
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Figure 5.1: RF Circuit Block Diagram. 


A. ANTENNA ELEMENTS 


The open-ended waveguide antennas that were used in the OSNS design are too 
wide to be used in the RSNS array. Since the element spacing for the RSNS array is 
small, micro-strip antennas were developed to replace the open-ended waveguide 
antennas. The stripline antennas are half-wave dipoles designed to resonate at 8.0 GHz. 
A schematic diagram of a micro-strip dipole antenna is shown in Figure 5.2. The width 
of each element is only 60 mils. The elements are placed into a linear array with the 
spacings set at 0.80" between the reference and modulus 17 elements and 1.70" between 
the reference and the modulus 8 elements as shown in Figure 5.3. These spacings were 


chosen over the original design spacings of 0.70" and 1.48", since the folding waveform 
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becomes distorted at angles greater than 60° from broadside. By compressing the RSNS 
dynamic range into the reduced field of view a higher angular resolution is obtained. An 
additional benefit of the increased spacing is that any possible mutual coupling effect is 
decreased. Antenna pattern measurements were taken of this array one element at a time; 


both phase and power were recorded. Figures 5.4 and 5.5 show the results respectively. 


DIPOLE ARMS 
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Figure 5.2: Micro-strip Dipole Antenna Element [from 5}. 
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Figure 5.3: Stripline Antenna Array [from 5]. 
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Figure 5.4: Phase Difference Between Modulus and Reference Signals. 
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Figure 5.5: Power Measured at the Mixer Inputs. 





B. ISOLATORS 


M2 Global S70124 isolators were used to further reduce mutual coupling effects. 
They are placed between the antenna elements and the first amplification stage. The 


specifications for the isolators are shown in Table 5.1. 





Frequency | Isolation ee eR VSWR Temperature Connector | Package 
Range =| (dB) Minimum Loss Maximum Range (°C) | Available Size (in) 
(GHz) (dB) Maximum 

| ee ee ee 








mO- 12.4 20 0.5 Zo -20 to +60 SMA 5 xX .90 xX .50 


Table 5.1: Isolator Specifications. 


After testing the isolators on the network analyzer they were found to meet all the 


specifications that were supplied. 
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c. AMPLIFIERS 


The system was first designed with only one stage of amplifiers. The one 
amplifier design was introduced to try and eliminate sources of phase errors that were 
present in the first prototype. The DBS DWT 18636 amplifiers were used. The goal of 
the amplification stage is to gain a constant power output to send to the mixers since any 
mismatch in power at the mixer inputs will produce phase errors. These microwave 
amplifiers have two operating regions: linear and saturation. The preferred region for 
this system is the saturated region. By operating in saturation, the output power of the 
amplifiers is independent of the input power. Figure 5.6 shows the power curves for the 


DBS amplifiers. 
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Figure 5.6: Power Curves for DBS Amplifiers. 


In order for the system to work well the input power needs to be at least -15 dBm. The 
initial transmitter distance was 19 feet. The input power for the amplifiers at this distance 
was not enough to drive the amplifiers into saturation. A stand was constructed to move 
the transmitter closer to the antenna elements. However, even with the transmitter at 6' 
from the elements, the power input to the amplifiers was only -17 dBm. This level was 
barely enough to saturate the amplifiers, so a two-stage amplifier design was utilized. 

The DBS amplifiers are used in the second stage of amplification because the 


power output levels are well matched and the three were phase matched to within 10°. 
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The first stage amplifiers were chosen based on their power curves. The two primary 
considerations were Maximum gain and how well the power curves matched. Figure 5.7 


shows the power curves for the three Avantek amplifiers that were chosen. 
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Figure 5.7: Power Curves for Avantek Amplifiers. 


Phase matching for this amplification stage is not as important as the second 
stage. The purpose for the first stage is to drive the second stage amplifiers into 
saturation. The gain for these amplifiers is approximately 40 dB in the linear region, 
which 1s more than enough to saturate the second stage amplifiers at a received power of 
-50 dBm. The overall effect of the two stages is that the signal output from the second 
amplification stage remains a constant regardless of the received power by the antenna 


elements. 


D. FILTERS 


The amplifiers used in this system produce harmonics. The harmonic levels are 
20 dB lower than the fundamental, but introduce phase errors in the system in excess of 
20°. In order to eliminate these errors, filters are used after each stage of amplification. 
AIRTRON band-pass filters, with a pass band of 4-8 GHz, are used after the first stage 
amplifiers and MICROLAB/FXR LA-9ON low-pass filters are used after the second stage 


4] 





amplifiers. Each filter has an attenuation of 55 dB or greater in the stopband, which 


effectively eliminates all signals except the fundamental. 


De MIXERS 


A mixer iS a common microwave component that is used to convert a higher RF 
into an intermediate frequency (IF). Typically, the RF signal is combined with a second 
signal, referred to as the local oscillator (LO), at a non-linear diode. The resulting IF is 
approximated by Wis = Oy - Wio [10]. The + + Gj. and other product terms that are 
generated in the mixer are removed by filtering. For the RSNS antenna, the mixers are 
used as phase detectors. The RF and LO inputs are the same frequency, therefore, the IF 


output is a DC value. 


IF=A, cos|(w,, — Ors t ~ (6,, — Oi )|+ A, cos|(w,, + Ov ) a (,, — Ors ) 


; (5.1) 
+ [higher frequency terms]. 
Since @,, = @,,, 
IF = A, cos(o,, — 0, )+ A, cos(2mt)+ (higher frequency terms ). (5.2) 


With proper filtering the 2@t and higher frequency terms can be eliminated. The 


remaining DC term depends solely on the phase difference (o,, — 4, ), 
IF = A, cos(o,, —9,,) - (5.3) 
Two mixers were considered for this system, the Mini-Circuits ZMX-10G mixer 


and the Anaren 70540 balanced mixer. Table 5.3 lists the specifications for each mixer. 
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Table 5.2: Mixer Specifications. 






For this application, the important specification is the IF frequency range. By having the 
range set below 2 it effectively filters the output of the mixer. 
The next important parameter that needs to be examined is the phase response of 


the mixer. Figure 5.8 shows the phase response of both mixers when connected to the 


network analyzer. 
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Figure 5.8: Mixer Phase Response. 


By comparing the curves, the Anaren balanced mixer was chosen since the Anaren 
mixer's response has a more symmetrical shape about the 0° axis and has its maximum 
and minimum values at + 90°. Also noticeable is a DC offset of approximately -40 mV. 
This offset can be compensated in the level and bias circuits at the input of the signal 


processing circuit. One drawback that was discovered during testing is the Anaren mixer 
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has a low saturation level. Figure 5.9 shows the first set of folding waveforms with the 


power level to the inputs at 14 dBm. 
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Figure 5.9: Folding Waveforms from Saturated Mixers. 


The maximum and minimum regions of the folding waveforms appear to be broadened, 
indicating that the mixer is saturated. Attenuators were added at the mixer inputs until 


the regions regained their predicted shape. 


F. MISCELLANEOUS COMPONENTS 


The phase difference between the reference signal and the channel signal is the 
basis for AOA estimations. The phase difference at broadside defines the relationship of 
the folding waveforms. In order to set the phase angles in each of the channels, the 
lengths of the semi-rigid coaxial cable were precisely cut. Because predicting how the 
phase of the signal is affected by the amplifiers and other components is difficult, phase 
adjusters were added. The phase adjusters are variable length coaxial cables that contain 
no dielectric medium. At 8.0 GHz, one turn of the barrel will change the phase of the 
signal by 5° giving the phase adjuster a range of approximately 60°. 

Attenuators are also added to each signal line. The purpose of the first set of 


attenuators is to balance the power between the reference signal and the modulus signal at 
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the output of the second amplification stage. A second set of attenuators is added to 
reduce the power input to the mixers, since the 20 dBm output of the second stage 


amplifiers saturates the mixers, causing phase errors. 


G. OVERALL DESIGN 


The system is assembled on a 1/4" brass plate measuring 13 in. by 13in. The 
plate serves two purposes. The first is structural: by attaching the components to the 
plate any vibrational errors that might be introduced are eliminated. The second function 
of the brass plate is as a heat sink. The amplifiers will operate at a high temperature, 
which may introduce phase errors into the system. The brass plate keeps the entire 
system at a constant temperature that is only slightly above room temperature, thus 
minimizing errors. All coaxial cables are cut to length. These custom fit lines minimize 
bends in the cables and the number of connectors that are needed to assemble the system. 


Figures 5.10 through 5.12 show the front, bottom, and top of the antenna system. 
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Figure 5.10: Front View of RSNS Antenna. 
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Figure 5.11: Bottom View of RSNS Antenna. 
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Figure 5.12: Top View of RSNS Antenna. 
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VI. RSNS SIGNAL PROCESSING CIRCUIT 


The second half of the DF system is the RSNS signal processing circuit. This 
circuit is responsible for converting the folding waveforms into angles of arrival. The 
circuit consists of 4 basic parts: level and bias amplifiers, comparator network, RSNS-to- 
binary converters and the timing circuit. Figure 6.1 is a block diagram of the entire 


circuit. 
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Figure 6.1: Block Diagram of RSNS Processor. 


The circuit gets its input from the balanced mixers and the folding waveform is 
sent through the level and bias amplifiers. This part of the system converts the folding 
waveforms (phase response from each channel) from a +320 millivolt signal to a 1 to 13 
volt signal. The signal is then sent to a comparator network. The comparator network 
operates like a thermometer, the higher the voltage of the signal, the more comparators 


that are turned on. The next portion of the system to receive the signal is the timing 
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circuit. This circuit is responsible for latching the comparator outputs to the EEPROMS 
to perform the RSNS-to-binary conversion. The conversion is simply a matter of 
comparing the number of comparators that are turned on in each of the channels with the 
RSNS code, generating a binary number between 0 and 63 that corresponds to an angle of 
arrival. 

The improvements made to the first prototype signal processing design include 
the addition of a timing circuit to latch the signal inputs to the EEPROMs (RSNS-to- 
binary converter) and constructing the preprocessing circuit on a printed circuit board 


(PCB). These modifications were added to increase the reliability of the system. 


A. LEVEL AND BIAS AMPLIFIER 


In order for the comparator network to function efficiently the + 320 millivolt 
signal needs to be conditioned. The first item to consider is that the signal takes on both 
positive and negative values. This changing polarity of the signal requires that the 
comparators have multiple power sources. To avoid multiple power sources a 7.0 volt 
bias is introduced. The second item to consider is that the range of the signal 1s only 640 
millivolts. For this case, the comparator voltage references would require very tight 
tolerances and very accurate comparators. This problem is alleviated by amplifying the 
signal so that its range is now 12 volts centered about a 7 volt bias. 

A weighted summing amplifier with an inverter is used to condition the signal. 
Figure 6.2 is a diagram of the weighted summing amplifier. The amplifier is a National 
Instruments LM-741CN Operational Amplifier. This amplifier has two inputs: the mixer 
output (Mod Signal) and V2 that acts as a bias voltage. R2 determines the amount of bias 
V2 provides. R, determines the influence generated by the mixer input. R3 sets the gain. 
The equation that governs the amplifier operation is the summing amplifier gain 


equation: 
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Nee Ns 
ve -a[ Bee] (6.1) 


With an +18 volt source, V2 is 18 volts and R3 is set to 1OkQ. The two remaining 
resistance values are determined by matching the mixer output maximum and minimum 
values to the desired outputs of 1-13 volts. The negative sign on the nght side of 
Equation (6.1) indicates that this is an inverting amplifier. The mixer output has a DC 
bias of -40 millivolts; however, when a 7 volt bias is introduced it becomes insignificant 
and is ignored. A phase detector output of zero maps to the middle of the shift and bias 


band: -7 volts. When V,=0, (6.1) reduces to 


V V 
R,=-R, <7 = ~10KSQ- = 25.722, (6.2) 


out 





The value of R, is determined from the maximum value of the mixer output. For the 
modulus 8 and 17 channels, the maximum mixer output 1s measured as 320 millivolts and 


the output is —13 volts. R, is determined as follows: 


Sede 0200 05207 
a a 


Vou, V2) (=13V,_18V ~ = 600HA 
R, R,} (10kQ  25.7k2 


Here R; is al kQ variable resistor set to 533 Q and R> is a 24 kQ resistor in series with 


= 533Q. (6.3) 





a1.6kQ resistor. 
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Figure 6.2: Weighted Summing Amplifier [from 5]. 


The output of the weighted summing amplifier must be inverted before 
transmission to the comparators. The inverter shown in Figure 6.3 performs this. No 
amplification is required in the inverting stage because the gain of the summing amplifier 


is correct . The equation that governs the output of an inverting amplifier 1s 


Ve ye (6.4) 





To obtain unity gain with inversion, R}=R2, 20k Q2 resistors are used. The output of the 
inverter is +1 to +13 volts as desired. This signal is provided to the comparator network 


to determine the level of each signal [5]. 
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Figure 6. 3: Inverting Amplifier [from 5]. 


DB: COMPARATOR NETWORK 


The comparator networks have 8 and 17 comparators in parallel. The circuit for 


each comparator is shown in Figure 6.4. 
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Comparator Circuits 
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Figure 6.4: Comparator Schematic and Pin-out Diagram. 


The shift and bias circuit output is the input to the non-inverting terminal. The second 
input is a reference voltage generated in the voltage divider network shown in Figure 6.5. 
Since the 18-volt power supply is available as the power supply for the comparator, it is 
used as the source for the voltage divider. The thresholds for the comparators, as 
explained in the previous section, range from 1 to 13 volts, so the conversion from the 


normalized thresholds to the actual thresholds is V,., =6V,.,, + 7. The normalized and 


actual threshold values are shown in columns two and three of Table 6.1. Applying 
Ohm’s law between the source and inverting terminal of the comparator, Equation 6.4 


shows the relationship between I, and Rg: 


= [ovis Yee (6.4) 


I, 
2 R, 
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The comparators (National Semiconductor LM311N) have a maximum input bias and 
offset current of 100 nA. This current is I, in Figure 6.5. In order to limit the input offset 
and bias currents to less than 1% of the current in the voltage divider, the total resistance 
of the network must not exceed 1.5M&2. However, the resistance values must be large 
enough such that the current required by the system remains small. For these reasons, Ry 
is chosen to be 100k 82. The total current for the voltage dividers is 27.5mA as shown in 
the lower left corner of Table 6.1. Kirchhoff’s Current Law applied at the node above R> 


yields Pt ee The only variables that remain are Rs and Re. Ohm’s law across these 


resistors dictates: 





V 
R,+R,=—™. (6.5) 


As an initial value, Rs is chosen to be 5% of Rg. Once initial values are determined, 
standard resistor values are chosen and Rs is made a variable resistor. The voltage 
divider sub-circuit is complete and provides a precise threshold voltage to the inverting 


terminal of the comparator [5]. 


a5 


Simulation Thresholds Gwen R4=100k and I2=100 nA 
Mod 8 Channel Resistance Values 


Normalized Comparator] R&R Use Ver Re RS Use Resistor 


1.11528 1.000E+05/ 812 314.723 


0.67370 
0.79802 
0.89516 
0.96183 
0.99573 


2.01116 
3 66655 
5.82943 
8.17052 
10.33341 
11.98881 
12.68471 


1.02559 
1.22905 
1.62902 


221190 
2 oaloe 
3.84141 

4.83255 
5.89750 

7.00000 
8.10250 
9.16745 
10.15859 
11.04217 
11.78810 
12.37098 
12.77095 
12.97441 


822 5.994E+02 
832 1.219E+03 
842 2.283E+03 
852 3'962E+03 
862 6.427E+03 
872 9.513E+03 
882 1.202E+04 


1 QOOE +05 
1.000E +05 
1.000E +05 
1.000E +05 


2.879E +02 
3.492E +02 
4741E+02 
6.676E +02 
3.370E+02 
1.293e+03 
1.749E +03 
2.322613 
3.033E+03 
3.902E+03 
4 948E +03 
6.177E+03 
7 568E +03 
9.051E+03 
1.048E +04 
1.165E +04 


812 
822 
832 
842 
852 
862 
872 


1.00E-07 
1.00E-07 
1.00€-07 
1.00€-07 
1.00E-07 
1.00E-07 
1.00E-07 
1.00E-07 


1.697E-04 
1.677E-04 
1.637E-04 
1.579E-04 
1.504E-04 
1.416E-04 
1.317E-04 
1.210E-04 
1.100E-04 
9.898E-05 
§.833E-05 
7.841E-05 
6.958E-05 
6.212E-05 
5.629E-05 
5.229E-05 
5.026E-05 
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1.687E-04 
1.598E-04 
1.432E-04 
1.216E-04 
9.819E-05 
7.657E-05 
6.001E-05 
§.105E-05 


1.696E-04 
1.676E-04 
1.636E-04 
1.578E-04 
1.503E-04 
1.415E-04 
1.316E-04 
1.209E-04 
1.099E-04 
9.888E-05 
8.623E-05 
7.831E-05 
6.948E-05 
6.202E-05 
5.619E-05 
5.219E-05 
5.016E-05 





6294 455 

1.199E +04 
2.438E+04 
4 565E +04 
79245.061 
1.285E +05 
1.903E+05 
2.404E +05 





A=10 
B=11 
C=12 
D=13 
E=14 
F=15 
G=16 
H=17 


Table 6.1: Voltage Divider Calculations. 


+18 V 





Figure 6.5: Voltage Divider Circuit [from 5]. 


To improve the comparator's performance a Schmitt tngger is added. R7 is 
introduced to add hysteresis to the comparator characteristics. This gives the comparator 
two threshold values: Vr, and Vry. Figure 6.6 shows the transfer characteristic of the 
comparator circuit. Vitsesh 18 set by the voltage divider and L is the comparator output 


signal. Vy, and V7y can be found by equations 


ie. 
e R,(R, +R.) (6.6) 


— 
TL thresh + R, 


and 


thresh 


his 
_y_ | ——___7RalRs + Ro) | (6.7) 
R, 


The power supply of the LM311 comparator sets the output signal, L=Vcc, in this case 
the L,= 18 volts and L. = -18 volts. 
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| 
| | 
, Hysteresis 


| 
! 
Figure 6.6: Comparator Characteristic with Hysteresis. 


A small capacitor, C, is added across the feedback resistor to ensure fast transitions and 
to prevent multiple pulses as the signal passes though the threshold voltages. The value 


of 10 M2 is selected for Ry and 100 pF for C [11, 12]. Table 6.2 lists the threshold 


values for each comparator. 
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Table 6.2: Comparator Threshold Settings. 


When the shift and bias signal is larger than Vry, the comparator is ON and 
remains on until the signal falls below Vy,. In ON state, no current is accepted at the 
Output, so the 5 volt signal above the output resistor has no voltage drop. The 5 volts is 
impressed on the output. This voltage is passed to the input of the latch. If the reference 
voltage exceeds the shift and bias signal, the comparator is OFF. In this state, the output 
sinks current, and the 5 volts are dropped across the output resistor. The output of the 


comparator is shorted to ground through the opamp and passed to the input of the latch. 


OF LATCHING THE COMPARATOR OUTPUTS 


A D flip-flop is added to the prototype design to latch the data to the EEPROMs. 
A set of latches is placed on the inputs for each of the two EEPROMs. A total of 39 
latches are required for the design. Six 74HC273N octal D type flip-flops are used as 
latches. Figure 6.7 shows the pinout for the 74HC273 chips. 
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Figure 6.7: 74HC273N Octal D-type Flip-flop. 


The V,, and CLR pins are set to 5 volts, D pins are the inputs, Q pins are the outputs and 
the CLK pin receives a signal from a 555 astable multivibrator. This pinout simply stores 
the input value on the D pin for one clock cycle. By latching the value for one clock 
cycle the inputs to the EEPROMs are stabilized. 

The clock signal is generated by 555 timer circuit shown in Figure 6.8. Ra, Rg, 


and C control the period and duty cycle. The design equations governing this circuit are 


Period = 0.69C(R, +2R,) (6.8) 
and 
R, +R 
Duty Cycle = pow allie (6.9) 
R, +2R, 


To select a clock signal of 1 kHz and a duty cycle of 50%, Ra, Rg, and C are set to 1.3 
kQ2, 6.2 kQ, and 100 nF respectively [12]. 
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Figure 6.8: 555 Timer Circuit. 


D. RSNS-TO-BINARY CONVERTER 


The 25 comparator states must be processed into a 6-bit angle of arrival. The 
EEPROM performs this mapping function. An ATMEL Model AT28C256 is used. The 
EEPROM is a small logic block that converts 15 input bits to an 8-bit output via its 
internal programming. Because of the 15 input limit, two EEPROM are used to process 
all 25 comparators. The “Thermometer Code” programmed into the EEPROM is shown 
in Table 6.3. The thermometer code is the string of instructions for the EEPROM that 
maps input to output. With two EEPROMs, the thermometer code has two stages as 


shown on the circuit diagrams of Figure 6.9 and Figure 6.10. 
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Figure 6.9: Schematic Diagram of Mod 17 Comparator Output to EEPROM#2 Input. 
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Figure 6.10: Schematic Diagram of Mod 8 Comparator to DIO Board Wiring. 
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The first stage encodes the 15 least significant (lowest threshold) comparators of 
the modulus 17 channel into a 4-bit representation (thermometer-to-binary conversion). 
The first column of Table 6.2 is the reported angle bin. The second column is the number 
of comparators active in the modulus 17 channel in each bin. The fourth column is a 
binary representation of the state of the input pins for the first EEPROM, and the fifth 
column is the hexadecimal representation of that binary number. The sixth column is the 
decimal output desired from EEPROM #1 for the given input; it is the decimal number of 
number of active comparators. The seventh column is the hexadecimal representation of 
that digital number. The result is that the number of active comparators is converted to a 
4-bit number. The shift and bias circuit and comparators for Modulus 17 and EEPROM 
#1 are physically located on the same circuit board. The output of EEPROM #1 and 
comparators 16 and 17 must be passed to EEPROM #2 on the second circuit board. 


Figure 6.9 illustrates how this signal 1s passed to EEPROM #2 [5]. 


E. CIRCUIT SIMULATION 


A schematic capture and simulation software package is used to test and validate 
the entire circuit design. Multisim ™ from Electronic Workbench was chosen because it 
is powerful, easy to learn and use. No knowledge of SPICE or other simulation language 
is needed to conduct basic simulations. The only knowledge that is necessary is basic 
circuitry design. The interface is graphical and very intuitive. 

Validation for each of the parts of the RSNS processor is done by selecting all the 
individual resistors, capacitors, integrated circuits, power supplies, and connecting the 
correct pins together. The program allows the user to add various measuring equipment 


to the circuit to monitor performance. Items such as oscilloscopes and voltmeters are 
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very useful to verify and tweak the initial design. The built-in function generator was 
used to provide the mixer output signal and each of the EEPROM inputs was monitored 
to ensure that the comparator network operated properly. 

A simulation of the complete system was not completed. The difficulty arose 
from the EEPROM. No standard SPICE model exists for EEPROMs, so a model needed 
to be created. Included in Multisim™ is a component creator that allows the user to 
generate their own unique components. This ability is limited to basic components such 
as diodes and opamps. In order to create a model for the EEPROMs a VHSIC Hardware 
Description Language (VHDL) add-on package is required. This language allows a 
programmer to model all levels of circuit design. It supports behavioral description of 
hardware from the digital system level to the gate level. The language is able to take a 
black-box approach to designing a system since all that is required in input/output 
information. This makes the language technology independent and extremely flexible. 

Three approaches were considered to simulate the behavior of the EEPROM. The 
first method is the easiest conceptually and most direct method. The truth table is entered 
into the VHDL code as a constant [13]. For smaller truth tables this method gives the 
programmer direct control of all input and output bits. For this system, there are 15 input 
bits that need to be accounted for thus requiring a truth table that has 2’° lines. The 
second method uses a programmable logic array (PLA) instead of truth tables. The table 
is replaced with an AND array that realizes selected product terms of the input variables 
and an OR array that operates on the product terms to form an output. This approach 
eliminates the need to include every possible input combination, thus reducing the 
number of lines of code. The drawback to this method though is that equations need to 
be developed and entered into Karnaugh maps, and a separate VHDL library file needs to 
be created to use this method. The third approach is the Case Model [14]. This approach 
is usually used to realize a multiplexer, however, it can be manipulated to simulate the 
operation of a ROM. The case statement allows the programmer to directly read the 
input bits and convert them into an output. This is similar to programming in a truth table 
except that only the required lines from the table are included instead of all the lines. 


Appendix A is the code developed to simulate the EEPROM. 
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Once the VHDL code was developed, it can be imported to Multisim™ and 
attached to a component. The component can then be added to the circuit for simulation. 
When the EEPROM components were added to the simulation, however, the Multisim™ 
package crashed. After troubleshooting, it was discovered that the VHDL interface for 
Multisim™ was not functioning. Electronics Workbench is developing a newer version 
of their design software to correct this bug, however, the fix was not available in time for 


the submission of this work. 


F. PRINTED CIRCUIT BOARD CONSTRUCTION 


In order to produce a PCB for the system, the entire circuit needs to be 
constructed in a program capable of generating a Gerber file (drawing file used by PCB 
milling machines). The Multisim/Ultiboard/Ultiroute™ software package produced by 
Electronics Workbench was selected to do this design. Multisim™ is a schematic capture 
and simulation program. It 1s used to design and test circuits. Once a design is complete, 
the program will generate a parts list and netlist, which is necessary to lay out the PCB. 
The Ultiboard program converts the parts list into a nearly complete scale drawing of all 
the components and the board itself. A board size is selected and the components are 
arranged to fit. Once the component layout is complete, the Ultiroute program converts 
the netlist into a network of traces. It does this conversion automatically and optimizes 
the design to reduce the number of traces required in the system. For complicated boards 
the routing process will invariably miss a couple of pins. This is corrected by manually 
adding the traces. 

The next step in the process 1s to generate a Gerber file for production. There are 
a few different types Gerber files that Ultiboard™ will generate. The equipment that is 
used to manufacture the PCB dictates the type that needs to be generated. Once the 
Gerber file is created, it is sent to a manufacturing site. The PCBs that were constructed 
for this system were milled by the NPS Physics department. The department owns a 


small milling machine capable of manufacturing 11 in. by 9 in. circuit boards. Because 
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of the size limitations of the milling machine, two circuit boards were designed. Once 
the boards were milled, ground and continuity checks were completed. When each of the 
boards was certified, they were stuffed and soldered. Tables 6.4 and 6.5 list all the 


materials required to assemble the printed circuit boards. 


Bill of Materials Modulus 8 PCB 















Quantit Description Reference ID Package 
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Table 6.4: Bill of Materials for Modulus 8 PCB. 
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Bill of Matenals Modulus 17 PCB 


Quantit Descniption Reference !D Package 
[SsSAPRESISTOR, 6.2kohm Ricca RESUS ee 
| _AJRESISTOR, 1.3kohm RiGee ee REGIIG. ee 


17/RESISTOR, 10Mohm 






R116, RI15, R114. 





R113 
pC IICONNECTORS, HDR2X7 TT CONIB DRO 
[TIMER UMSS5CN NE 
39|CAPACITOR, 100nF cap4 
, C5, C6, C35, C36, 
. C37, C38, C34, C33 
, C32, C31, C38, C2 
3, C28, C27, C26, C 
25, C24, C23, C22, 
C21,,C20, C19, C16, 
C17, C16, C15, C14 
, C13, C12, C11, C1 
0,C3, C8, C7 
PS SAU S RANG OTON a NO 
17|CAPACITOR, 100pF C75, C74, C73, C72, Cap2 
RS Ss S77 ae Sara ee 


17| COMPARATOR, LM311N 


RESISTOR, 100kohm 





S 





RESISTOR, 22kohm RESO.5 
R68, R654, R60, R55 
_ R51, R54, R43, R4 


7, R37, R34, R21,R 






Dic re olo 
mi = m aa 
Bio ¥ ¢o 
in| O us 6 
= 

AL!) 


25, R28, RIO 
a) SiGe i res ae EE eae 
ee 
R24, Reg 
SS) (Sonn SS iS aa SS ee ae 
| 2IRESISTOR, 20kohm RAR RESO 
| IIRESISTOR, 10kohm RESO 
| IRESISTOR, 24kohm IR RESOS 
J CIRESISTOR, 76ohm RG RESO 
| AIPOTENTIOMETER, 100K UN RZ RD, RES, RET LINPOT 
eee RESISTOR, 200kohim 55 RE 
| 2IRESISTOR, 150kohm RRS RESO 
Wi Roa a eee eS a eee | 
J IRESISTOR, 120kohm RRS 
| IRESISTOR, 75kohm RAB RESO 
See a SS ara ae 
RD 0S 7, - (== 
Re SCE i - aee eae 
| TIRESISTOR, 12kohm RR RES 
| AJRESISTOR, 18kohm RG RESOS 


Table 6.5: Bill of Materials for the Modulus 17 PCB. 
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After the board is stuffed and soldered, power is connected. The component 
connections are tested and the trim potentiometers are set to provide the correct bias and 
amplification to the signal and provide the proper reference voltages for the comparators. 
Now the boards are ready for operation. The Modulus 17 and 8 PCBs are shown in 
Figures 6.11 and 6.12 respectively. Appendix C gives a detailed narrative of how to 


manufacture a PCB. 
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Figure 6.11: Modulus 17 Printed Circuit Board. 
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Figure 6.12: Modulus 8 Printed Circuit Board. 
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Vil. TEST AND EVALUATION OF PROTOTYPE HARDWARE 


A. TEST PROCEDURES 


The prototype antenna was tested in the Naval Postgraduate School Microwave 
Anechoic Chamber. The major components and communication paths used for testing 
are shown in Figure 7.1. The test system architecture centerpiece is a Pentium I 
computer. One of two LABVIEW programs controls the test via a General Purpose 
Interface Bus (GPIB), the sertal port, and a PC-DIO-96 board. The signal generator that 
produces the radiated signal is located on the same rack as the Hewlett Packard 8510 
Vector Network Analyzer (VNA). In the anechoic chamber, a standard gain feed hom 
transmits the microwave signal generated by the signal generator. The antenna is placed 
on a pedestal (Figure 7.1) approximately 19 feet from the transmit hom with the 
microwave circuit mounted on the back of the array ground plane. A servomotor is 
controlled through the computer serial port and rotates the pedestal, which is capable of 
steps as small as 0.1°. The microwave circuit continuously analyzes the received signal, 
and the mixer output 1s sampled at each resolution step. The mixer output is carried out 
of the anechoic chamber on RG-59 coaxial cables and either connected to the RSNS 
processor or two HP 3478A multimeters. The 6-bit output of the RSNS processor is 
transmitted to the computer via the DIO board. Two HP 3478A digital multimeters 
measure either the mixer output or the shift and bias circuit output. The computer 


samples the multimeters through the GPIB bus. 
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Figure 7.1: Anechoic Chamber Setup. 


The data collection and processing technique is illustrated in Figure 7.2. Two 
types of tests are conducted depending upon the parameter of interest. The first test 
measures the output of the mixer, and other test measures the output of the shift and bias 
amplifier and the binary output of the RSNS processor simultaneously. One LABVIEW 
program is written for each test. When the mixer output is measured, the procedure 
follows the left path on Figure 7.2. The LABVIEW VI code DF_mm_R2 samples the 
multimeters at each servomotor step. The data is saved in column format. The four 
columns of data are (1) antenna position, (2) channel | (modulus 8) multimeter output, 
(3) antenna position, and (4) channel 2 (modulus 17) multimeter output. The data is 


analyzed with Microsoft Excel using the “text-tab delimited” format. Manual 
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normalization of the data is performed. The minimum and maximum values of the phase 
detector folding waveforms are located. A column is inserted between the second and 


third data columns. The normalized phase detector output Ni norm 1S Calculated using the 


equations: 
ie Oe tty) ea (7.1) 
AN = (Nnax Se ee (7.2) 
n-7 
| es TS 
Pom An ( ) 


The normalized output for each channel is placed in the new third column (modulus 8) 
and the sixth columns (modulus 17). The normalized phase detector output is saved in 
tab delimited text format with a “.dat” extension. This is the standard input format for 
the MATLAB processing programs provided in Appendix A. The measured phase 
(mixer Output) is visually compared to the simulated data using the rsns&.m and rsns17.m 
programs. Phase_error.m calculates and plots the phase error between a simulated and 
normalized measured mixer output. The normalized mixer output is passed through a 
simulation of the digital processing in Bin_Pop.m. 

The output of the digital circuit follows the path on the nght of Figure 7.2. The 
shift and bias circuit output 1s connected to the multimeters. The binary output of the 
digital circuit is connected to the computer via the DIO board. The LABVIEW VI code 
DF5_DIO_R2 samples the multimeters and DIO board at each servomotor step. The data 
is again saved in column format. No normalization is required for this data, but it must 
be saved in tab delimited text format with a “.dat” extension. Analysis of the data is 
performed by one program, AOA_map.m. This program plots the folding waveforms and 
the bin mapping. It converts the bin numbers into an estimated angle of arrival and plots 


the transfer function. 
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Phase Shift and EEPROM #2 
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Figure 7.2: Anechoic Chamber Control and Recording Equipment. 
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B. RESULTS 


The RSNS direction finding system was built and tested. The array design is 
scaled so that the field of view is 60°. The tests were conducted at 8 GHz at a resolution 
of 0.1°. Data from the mixer outputs, RSNS processor, and level and bias amplifiers 
were recorded. 

The folding waveforms for both moduli are shown in Figures 7.3 and 7.4. The 
blue curves are the ideal waveform and the red curves are the measured waveforms after 
they have been normalized. The green curve is the residual error of the measured signal. 
BinPop.m produces the simulated RSNS antenna transfer function using measured mixer 
outputs. Figure 7.5 shows this transfer function. Once it was shown that the transfer 
function was mostly continuous, a complete test of the system was conducted. Figure 7.6 
is the folding waveforms after the shift and bias amplifiers have processed them. The 
output of the RSNS processor is recorded through the DIO board and the AOA_map.m 


file produces the angle of arrival estimation of the entire system, shown in Figure 7.7. 
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Figure 7.3: Modulus 8 Normalized Mixer Output. 
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Figure 7.4: Modulus 17 Normalized Mixer Output. 
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Simulated RSNS Antenna Transfer Function Using Measured Mixer Outputs. 
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Figure 7.6: Level and Bias Amplifiers Outputs. 
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Figure 7.7: RSNS Antenna Transfer Function. 
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VII CONCLUSIONS 


A. CONCLUSIONS 


This research has introduced the Robust Symmetrical Number System and shown 
that it can be used to resolve ambiguities in phase sampling interferometry. A direction 
finding antenna architecture based on the RSNS has many inherent advantages over 
conventional direction finding techniques including: simple microwave beamforming 
network that incorporates wideband components, wide instantaneous field of view, and 
high resolution angle of arrival with a vanety of element arrangements having very short 
baselines. The design equations for comparator thresholds, broadside phase difference, 
and a scale factor are derived. The simulated transfer function demonstrates that the 
RSNS antenna will have a wide instantaneous field of view and high resolution. 

A second prototype array was built and its performance measured. Several 
improvements were made to the original RSNS prototype antenna. By adopting the 
scaled antenna carriage, adding isolators, and completing the ground plane with copper 
tape between the antenna elements mutual coupling effects are reduced. Mounting the 
microwave components on a brass plate has reduced errors contributed by vibrations and 
temperature. Tailor cutting all semi-ngid coaxial lines has reduced the number of 
connectors required to assemble the microwave circuit and has subsequently reduced the 
errors inherent in connectors. Matching amplifiers by amplification stage rather than by 
signal line has reduce relative phase errors between channels as well as better matched 
the power outputs of both amplification stages. Two printed circuit boards were designed 
and built for the RSNS processor. The printed circuit boards provide a more stable 
platform over the original design, which was assembled on breadboards. 

The measured transfer function of the new prototype exhibits the basic features of 
the ideal transfer function. The measured transfer function contains small encoding 
errors. These errors are attributed to phase errors in the microwave circuit and its 


components as well as threshold errors inherent in the RSNS processor. Comparing the 
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transfer functions of the first and second prototypes and the ideal case, the improvements 
in the system are revealed. Figure 8.1 shows the simulated transfer function of the first 
prototype RSNS antenna using the measured mixer outputs. Note that the large phase 
errors present (see Figure 4.7 and 4.8) cause a significant degradation in the antenna 


performance. 


Reported Angle of Arrival (degrees) 
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Figure 8.1: Simulated Transfer Function of the RSNS Antenna Using Measured Mixer 
Outputs (First Prototype) [from 5]. 


Figure 8.2 shows the simulated transfer function of the second prototype RSNS antenna 
using the improved mixer output response (see Figure 7.3 and 7.4). Note that this 
response matches the ideal response (shown in Figure 8.3) much more closely with only a 


few encoding errors present. 
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Figure 8.2: Measured Mixer Output Through Simulated RSNS Processor, Second 
Prototype. 
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Figure 8.3: Ideal Transfer Function. 
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The initial results from the RSNS digital processor shown in Figure 8.4 were not 
expected. The EEPROM programming will output zeros if the inputs do not match a 
known set of inputs from Table 6.3. This occurs in the regions where the estimated angle 
of arrival is -90°. Each set of components (i.e. program on EEPROM, latches, 
comparator thresholds, and bias and level amplifier performance) was individually 
verified for proper function; however, when the whole system is tested it does not 
entirely work. Errors in the circuit design may be attributed to the comparators 
themselves. Data sets from different runs were compared against each other to see where 
the errors occurred. It was noticed that the errors sometimes occurred at a specific angle 
and sometimes not. Using the anechoic chamber and testing a single angle verified this. 
Even though the comparator threshold's accuracy is improved with the Schmitt trigger, 
errors may still be occurring. To test this, the Labview code needs to be modified so the 
comparator states, input signal, and threshold values can be monitored. If it is 
determined that the comparators are not accurate enough, a precision comparator such as 


the CMP02 or CMP04 could be substituted for the LM311. 
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Figure 8.4: Measured Antenna Transfer Function. 
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B. FUTURE WORK 


The complete design of the RSNS direction finding system is an continually 
evolving process. A number of areas for further research exist. These areas include 
topics such as frequency agility, increased resolution, and operation in a multi-signal 
environment. 

It is desired to design a new system capable of handling a wide range of 
frequencies. The angle of incidence of the signal determines the phase difference 
between the antenna elements. This phase difference is linearly dependant upon the 
wavelength and therefore the frequency of the incoming signal. A relationship between 
the angular resolution and the dynamic range of the RSNS can be derived. Research into 
how to control the frequency effects on the RSNS antenna array is needed. 

Another improvement to the design that can be looked into is increasing the 
resolution of the system. There are several methods that could be used to accomplish this 
task. Some ideas are to increase the dynamic range of the RSNS by increasing modulus 
number in each channel or by increasing the number of channels. Preliminary data 
suggests that phase errors are more prevalent in higher modulus channels therefore this 
might not be a desirable method to accomplish increased resolution. Work needs to be 
conducted to determine the optimal method to accomplish this task. 

Finally, in order to be able to field a system in an uncontrolled environment, a 
multi-signal capability needs to be included into the design. It maybe possible to select 
individual frequencies by adding an adaptive filter or frequency selector to the front end 


of the system. 


c, CONCLUDING REMARKS 


The focus of this thesis research has been to improve the design and 
performance of a previous RSNS DF array. The design equations and operating 


parameters have been presented. The overall performance of the microwave circuit has 
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been improved and the phase errors reduced to a level to where the signal processor con 
effectively operate. Although the results from the RSNS signal processor are 
disappointing it has been shown that the RSNS signal processing technique is a viable 


technique for future DF applications. 
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APPENDIX A. VHDL CODE 


This appendix presents the VHDL code used to simulate both EEPROMs for the 
system. Each EEPROM needs two modules to function correctly. The first module is the 
actual code that is used to simulate the component. The second module is a test bench 


that is used to verify the functionality of the code. 


--EEPROM #1 Module 
-- Note: comments beginning with WIZ should be left intact. 
-- Other comments are informational only. 


-- Multisim VHDL software version: 3.75a 


library ieee; 
use lieee.std_logic_1164.all; 


-- use leee.numeric_std.all; -- Note: uncomment this if you use 
-- IEEE standard signed or unsigned types. 
-- use ieee.std_logic_arith.all; -- Note: uncomment this if you use 


-- Synopsys signed or unsigned types. 
entity EEPROM_1 is 
generic (D:time :=10 ns); 


Dore 
AO: in Dit; 
Al® in bit: 
AZ an bie: 
BS2hible 
A4: in bit: 
AS: in BLE; 
AG: in bit; 
A)? Ure DPE: 
A8: in bit; 


A9: in bit; 

AlO: in bit; 

All: in Dit; 

Al2: in Dit; 

ALS * yin bie. 

Al4: in bit; 

TOOs Out Dit; 

TOL. Sout Dit: 

162: out bit; 

103: Cut. bit 

Nee 

end EEPROM_1; 
-- From VHDL DeSign, J.R. Armstrong, pg. 320 
-- Model for combinational logic using CASE statement 
architecture EEPROM_1A of EEPROM_1 is 


begin 
Pl: process (A14,A13,A12,A11,A10,A9,A8,A7,A6,A5,A4,A3,A2,A1,A0) 
begin 


case A14&A13&A12&A11&A1 0&A9DS&ABEATEABEASEASEAZEAZEALEAD 1S 

when "000000000000000" => I03 <= '0O' after D; I02 <= '0' after D; 
LO1l <= 'O' after D; I00 <= '0' after D; 

when "000000000000001" => I03 <= '0' after D; I02 <= '0' after D; 
IO01 <= '0' after D; I00 <= '1' after D; 

when "000000000000011" => I03 <= '0' after D; I02 <= '0' after D; 
I01 <= '1' after D; I00 <= '0Q' after D; 

when "OOCOO0Q0Q000C00111" => 103 <= "0 after D> I02 <= ‘O° after D; 
IO1l <= '1' after D; I00 <= '1' after D; 

when "000000000001111" => 103 <= '0' after D; I02 <= '1' after D; 
101 <= 'Q’ after D; I00 <= '0' after D; 
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wnen "QQ00COO0G00J 1111". => 103 -<= ‘'O* ,atter D: 102 <= “1 atter 
IO). <=-"@" sateer -D- TO0:' <= ‘1° atter bp: 


when "000000000111111" => 103 <= ‘'0' after D; 102 <= '1' after 


10]. <= "1" after D: 100 <= °0”* after D; 
when "000000001111111" => I03 <= '0' after D; I02 <= '1' after 
TO1 “<= "1" after D2: 100 <= "1" after DD: 


when "000000011111111" => I03 <= '1' after D; I02 <= '0' after 


1Ol <=) "OO" atter DD 100 <=. "0" attrer Dp: 
when "“QOO000CGlI1II1111" => 103. <= "2" atter D: 102 <=. 'O'~atter 
101. <=" 0" -atiter D; 100 <=." 1" atter Dp; 


when "OOGOU001TU1ILI1TLI 1 =s 163 <=." after D: 102 <=) 0" “atter 


I01l <= ‘'l' after D; I0O0 <= '0' after D; 

when “OOO O11 TIa1 i" -=s-1O3-<= "1° atter D; 10Z i<=-"0" atte: 
IO] <= ‘1* after D; 100 <= ‘1 after D; 

wnen "OOOLTT1IttIiiiY => 163 <= "1" atter D; 102 <= *1" after 
101, <= “0" atter D; 100 <= "0" safter D; 

when "OOa Titres ss Os <= "1" -afttéeér D:/102 <=.*)" after 
IOl <= *0° atter DD: [100 <= *1'* after D; 

when "OLITIIII111i1111" => 103 <= *l after D; 102 <= "1" after 
TOL <=] ‘i -atter De 100 <=."0O° atter D; 

when "iii => 103 <=" 1.’ “atter De 102 <=)" 1" atcer 
10d <=" 1" areer DD; lO0-<- 1 ateer: 1); 

end case; 
end process Pl; 
end EEPROM_1A; 


mm re ar wwe we we we we we we we ww we we ee we ew a eae we i a ee 


--Test Bench Module for: EEPROM_1 
-- Note: comments beginning with WIZ should be left intact. 
~~ Other comments are informational only. 


-- Multisim VHDL software version: 3.75a 

library ieee; 

use ieee.std_logic_1164.all; 

-- use leee.numeric_std.all; -- Note: uncomment this if you use 


-- IEEE standard signed or unsigned types. 


-- use ieee.std_logic_arith.all; -- Note: uncomment/modify this if you use 
-- Synopsys signed or unsigned types. 

use std.textio.all; 

use work.all; 


entity TESTBNCH is 
end TESTBNCH; 


architecture stimulus of TESTBNCH 1s 
component EEPROM_1 is 
port. ( 

20.) in Bit; 
ALe In Bit? 
A2: in BIT; 
A3: in BIT; 
A4: in BIT; 
AS: in BIT; 
AG: in) BIT; 
A7: in BIT; 
A8: in BIT; 
A9: in BIT; 
ALO: in BIT: 
Ale in BIT 
AlZ: in BIT; 
ALS FBT s 
Ai4:\ in BIT; 
TOO: out BIT; 
LOl= out BIT; 
1602s" Out BL; 
103% Sut BIT 
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ie 


end component; 
constant PERIOD: time 


-- Top 
Signal 
signal 
Signal 
Signal 
Signal 
Signal 
signal 
Signal 
signal 
Signal 
signal 
Signal 
signal 
Signal 
Signal 
Signal 
Signal 
Signal 
Signal 


signal 


begin 


SO2ns - 


level signals go here... 


AO: 
Al: 
A2: 
A3: 
A4: 
A5: 
A6: 
A7: 
A&: 
Ag: 
Al0Q: 
All: 
Ra 
Al3:; 
Al4: 
IOO0: 
TOA: 
ra2. 
LO3: 


aH Bid Bae 
BL: 
BIT: 
BIT; 
po lel ey 
bo Bl aes 
Bit; 
Bi. 
Bit: 
BIT; 
BIT; 
BIT; 
BLT; 
BLT; 
BLT; 
Ble; 
Bit; 
BIT; 
BLL, 


Gone: boolean := 


false; 


DUT: EEPROM_1 port map ( 


STIMULUS1 : 


): 


begin 


AQ 
Al 
A2 
A3 
A4 
A5 
A6 
A7 
A8& 
AQ 


AO, 
Al, 
A2, 
A3, 
Ad, 
AS, 
A6é, 
A7, 
AS, 
AQ, 
Al0, 
All, 
Al2, 
Al3, 
Al4, 
Too, 
LO, 
IO2, 
I03 


f 


process 


Sequential stimulus goes here... 


Sample 


A A 
ooo 


“O4 
oa @ 
La Gib 
Li @ Th 
eo 8 
"Gis 
et 


A AAA 
nou ue ou 


A A 
nou 


A 
nou 


A 


Oa 
WiC 
a8 Bae 


Stimius =. < 


after 
after 
after 
after 
after 


aoaooo oO 


ns, 
ns, 
ns, 
ns, 
ns, 


vel 
Dale 
5 Re, 
Ue lg 
yo les 


after 
after 
after 
after 
after 


100 
200 
300 
400 
500 


ns; 
ns; 
nS; 
ns; 
ns; 
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AlGv<=."0"; 


All <a= "oO" 
Al2 <=.) 0"; 
Bis ee Os 
A1l4 <= 'OQ'; 
wait; -- Suspend simulation 


end process STIMULUS1; 


end stimulus; 


—— ecm ee me ee me me ce ce cm wm we we we ww we eee ee ea ee ei ei ei es 


--EEPROM #2 Module 
-- Note: comments beginning with WIZ should be left intact. 
-- Other comments are informational only. 


-- Multisim VHDL software version: 3.75a 


library ieee; 
use ieee.std_logic_1164.all; 


-- use ieee.numeric_std.all; -- Note: uncomment this if you use 
-- IEEE standard signed or unsigned types. 
-- use ieee.std_logic_arith.all; -- Note: uncomment this if you use 


-- Synopsys signed or unsigned types. 
entity EEPROM_2 is 
generic (D:time :=10 ns); 


port 4 
-- WIZ BEGINPORTS (Entity Wizard command) 
AO: in BIT; 
Al: in BIT; 
AZ = -in Bit: 
Asc. dnvBiT: 


A4: in BIT; 

AS: in BIT; 

Ab:> in. BIT; 

A?is in BIT; 

A8: in BIT; 

AS: in BIT; 

ALO2 in BET - 

All: in BIT; 

Al2Z. an BIT 

Al3: in BIT; 

LOO: out BIT: 

LOL? 7OUt- BIT; 

I02: out BIT; 

1O3:.-out. Birr: 

I04: out BIT; 

1O5:.<out- BIT 

-- WIZ ENDPORTS (Entity Wizard command) 
ie; 

end EEPROM_2; 
-- From VHDL Design, J.R. Armstrong, pg. 320 
-- Model for combinational logic using CASE statement 
architecture EEPROM_2A of EEPROM_2 is 


begin 
Pl: process (A13,A12,A11,A10,A9,A8,A7,A6,A5,A4,A3,A2,A1,A0) 
begin 


case A13&A12&A11&A10&A9&A8EA7TEABEASSASEAZEAZEALEAD IS 

when "00000000000000" => I05 <= '0' after D; I04 <= '0' after D; 
103 <= 'O* after BD; 102 <= "0" atter D; 101 <=," 0" after -D: 
Io0 <= '0' after D; --0 

when °00111111001010° => IOS <= '0’ after D; I104 <= '0' after D; 
I03 <= '0' after D; I02 <= '0' after D; I01 <= '0' after D; 
I0o0 <= '0' after D; --0 

when °00011111001010" => I05 <= '0’' after D; I04 <= '0Q' after D; 
103 <= '0" after D; 102 <= ‘0* after D; I01 <= ‘0° after D; 
I00 <= '1' after D; --1l 
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when *"00011111001011" => IO5 <= 
I03 <= '0Q' after D; I0O2 <= 
I00 <= '0' after D; --2 
when *"00001111001011" => IOS <= 
I03 <= '0Q' after D; IO2 <= 
IOO0 <= '1' after D; --3 
wnen *0000127211001100" => 305 <= 
I03 <= '0' after D; I02 <= 
IO0O0 <= '0' after D; --4 
when *00000111001100" => IOS <= 
IO3 <= '0' after D; I0O2 <= 
100 <= '1' after D; -=5 
when "00000111001101" => IOS <= 
I03 <= '0' after D; I02 <= 
IOO0 <= '0' after D; --6 
when "00000011001101" => IOS <= 
I03 <= '0' after D; I02 <= 
I00 <= '1' after D; --7 
when *"00000011001110" => IOS <= 
103 <= 'l" after D; I02 <= 
IO0 <= '0' after D; --8 
when "00000001001110" => I0O5 <= 
103 <='1l" atter D; 102 <= 
I0o0 <= '1' after D; --9 
when "00000001001111" => I05 <= 
I03 <= '1' after D; I02 <= 
I00 <= '0' after D; --10 
when "00000000001111" => IOS <= 
IO3 <= '1' after D; I02 <= 
IOO0 <= '1' after D; --i11 
when "00000000011111" => IOS <= 
IO3 <= '1l' after D; I02 <= 
Io0 <= 'Q' after D; --12 
when *"00000001011111" => IOS <= 
IO03 <= '1' after D; I02 <= 
I00 <= '1' after D; --13 
when "00000001111111" => IOS <= 
L03 <= '1* atter D; I02 <= 
IOo0O <= '0' after D; --14 
when "00000011111111" => IOS <= 
IO03 <= '1' after D; I02 <= 
IOO0 <= '1' after D; --15 
when *00000011011111" => IOS <= 
I03 <= '0' after D; I02 <= 
I00 <= '0' after D; --16 
when “OOGOQ1L1I1O1LI111" => I05 <= 
IO3 <= '0’° after D; I0O2 <= 
IOO <= '1' after D; --17 
when *"00000111001111" => I05 <= 
I03 <= '0' after D; I02 <= 
IOo0 <= '0O' after D; --18 
when *"00001111001111" => IOS <= 
IO3 <= '0O' after D; I02 <= 
I00 <= '1' after D; --19 
when *00001111001110" => I05 <= 
IO03 <= '0O' after D; I02 <= 
IOO <= '0' after D; --20 
when *00011111001110" => IOS <= 
I03 <= '0' after D; IO02 <= 
IOO <= '1' after D; --21 
when *%00011111001101" => IOS <= 
I03 <= '0' after D; I02 <= 
IOO <= '0' after D; --22 
when “OOLITI1I00LIOl" => 105 <= 
103 «= ‘0 after D;-102 -<= 
I0oOQO <= '1l' after D; --23 
when “OCLLIITIOOLIOO" => 105 <= 
103 <= '1' after D; I02 <= 
I00 <= '0O' after D; --24 
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after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


after D; 
after D; 


104 <= 
101 <= 


LO4. == 
Iol <= 


104 <= 
LO <= 


I04 <= 
I01 <= 


I04 <= 
Iol <= 


I04 <= 
TOL <= 


104 <= 
50 @ alee 


I04 <= 
101i <= 


I04 <= 
I0i <= 


14 <= 
IOl <= 


I04 <= 
1Ol <= 


I04 <= 
TOL <= 


L104 <= 
IOl1 <= 


I04 <= 
I01 <= 


104 <= 
I0l1 <= 


I04 <= 
Iol <= 


104 <= 
LOl> <= 


I04 <= 
I0l1 <= 


104 <= 
I0i1 <= 


I04 <= 
IOl1 <= 


164 -<= 
LOL <= 


164 <= 
Iol <= 


104 (<= 
Iol <= 
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after D; 
after D; 
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after D; 
after D; 
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after D; 
after D; 
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after D; 
after D; 


when "01111111001100" 


I03 <= '1' 
lOO <=." 


wher “Gli 7T111001011" 
after D; 102 <= 


TOS -<= 1 1;' 
1O0 <=" 0" 


when °11111111001011" 


LOS ncaa 
LOOe<=.V 4" 


when "11111111001010" 


Los <= 01? 
LOO. <=" "6* 


when "01111111001010" 


I03 <= '1' 
LOO <= "1" 


when *%01111111001001" 


103) <= 11h’ 
Io0 <= '0Q' 


when "00111111001001" 
areer D; 102 <= 


I03 <= '1' 
IoO0 <= '1' 


when "00111111001000" 
atter D; £O2 <= 


103 <= -*0* 
LOO = 20) 


when "00011111001000" 


163: <= "0° 
TOO <= 5-12 


when "00011111000111" 


EO3 “<=. 70" 
I0oO <= '0Q' 


when *%00001111000111”" 
after D; IO2 <= 


IOs <= "0" 
I0O0 <= '1' 


when *"00001111000110" 
after D; I02 <= 


103-<=~ 0° 
TOO" <= +*0* 


when "00000111000110" 
after D; IO2 <= 


LOS) <= 5-0" 
Io0 <= '‘'1' 


when "00000111000101" 


I0o3 <= '0' 
TOO. <S0° Oe 


when "00000011000101" 
after D; IO2 <= 


103) <=. 7.0" 
I0o0 <= ‘'1' 


when "00000011000100" 


LO3. <=! 1" 
I0oO0 <= 'Q' 


when *"00000001000100" 
after D; I0O2 <= 


1O3a <= 
100 <= * 1" 


when "00000001000011" 


103. <=::" 1’ 
100 <= "0" 


when "00000000000011" 


103 <= * 1’ 
POOt=<= ac! 


when *%00000000000010" 
after D; IO2 <= 


TO3¢ <=)" 
I0o0 <= '0' 


when "00000001000010" 


TOS sea ee 
I00 -<s* "1" 


when "00000001000001" 
after D; I02 <= 


TOs <=" 1! 
I0o0 <= '0Q' 


when "00000011000001" 


Os) <= ""1" 
I00 <= '1' 


I04 <= 
LOl <= 


I04 <= 


after D; IOl <= 


I04 <= 
TOl <= 


LO4 <= 
10] <= 


104 <= 
Iol <= 


after D; I04 <= 


Lod <= 


I04 <= 


after D; IOl <= 


I04 <= 
Iol <= 


I04 <= 
Iol <= 


I04 <= 


after D; IoOl <= 


I04 <= 
Iol <= 


I04 <= 


after D; IOl1 <= 


I04 <= 
P01 <= 


I04 <= 


after D; IOl <= 


I04 <= 
TOl. <= 


[O44 <= 
Iol <= 


after D; I04 <= 


LO <= 


I04 <= 


after D; IOl1 <= 


after D; I04 <= 
after D; IOl1 <= 


I04 <= 


after D; IOl <= 


after D; I04 <= 


Iol1 <= 


I04 <= 
Iol <= 


after D; I04 <= 


I0l <= 


‘l' after D; 
"O" ‘afeer D: 


‘l' after D; 
'l' after D; 


‘1* after D; 
‘1l' after D; 


ml atter Dp: 
'O' after D; 


"LY scatter Ds 
'O' after D; 


‘l' after D; 
‘1l* after D: 


‘1’ after D; 
"l* after D- 


'O' after D; 
"QO" “after D: 


'O' after D; 
"QO" after D:; 


'O' after D; 
“1” aftter D; 


'O' after D; 
"1" atter Db; 


7O* -akeer Ds 
‘0 “ealeer-p: 


'O' after D; 
'O' after D; 


'O* after D; 
'1l' after D; 


"0 satter “Dp; 
'l' after D; 


'O' after D; 
'O arter “D- 


"O° after D: 
nO) Sat cer <p: 


'O' after D; 
‘1' after D; 


'O' after D; 
*l* after D; 


'O' after D; 
'O' after D; 


"0" after D: 
"O° after D; 


'O' after D; 
"l* atter D; 


"0. atter Db: 
"l’ after D; 


when °00000011000000" => I05 <= '1' after D; I04 <= '1' after D; 
I03 <= '0O' after D; I102 <= ‘'0' after D; I0Ol1 <= '0' after D; 
I00 <= ‘0' after D; --48 

when °00000111000000" => I05 <= '1' after D; I04 <= '1' after D; 
103) <= ‘0° after DD; 102><= "0" after D> 10] <=.*0* after D; 
IO0 <= '1' after D; --49 

when "00000111000001" => I05 <= '1' after D; I04 <= '1' after D; 
L@35<=" 0" atter.D;-102(<= "0" atter Dp; 101-<=''1" atter D; 
IO0 <= '0' after D; --50 

when "00001111000001" => IOS <= '1' after D; I04 <= '1' after D; 
I03 <= '0' after D; I02 <= '0' after D; I01 <= ‘1' after D; 
IO00 <= ‘1' after D; --51 

when "00001111000010" => I0S5S <= ‘1' after D; I04 <= '1' after D; 
1O3 <= "0 watter Ds 102 <= 1 after D; 101 <= "0" aiter D:; 
I0o0 <= '0' after D; --52 

when °00011111000010" => I0S <= '1' after D; I04 <= '1' after D; 
1O3-<= "0" ~“atter D; 102 <= "1" atter D; 10] <= “0 -atter D; 
I00 <= '1' after D; --53 

when °00011111000011" => I05 <= '1' after D; I04 <= ‘1' after D; 
I03 <= '0O' after D; I02 <= ‘1' after D; IO1 <= ‘1' after D; 
IOO0 <= '0' after D; --54 

when "00111111000011" => I05 <= ‘1' after D; I04 <= ‘1' after D; 
I03 <= '0' after D; I02 <= ‘1' after D; IO1 <= ‘1' after D; 
I0o0 <= '1' after D; --55 

when "00111111000100" => I05 <= ‘'1' after D; 104 <= '1l' after D; 
I03 <= '1l' after D; I02 <= '0' after D; IO1 <= '0' after D; 
IoO0 <= '0' after D; --56 

when "01111111000100" => I05 <= ‘1' after D; I104 <= '1°* after D; 
I03 <= '1l' after D; I02 <= '0' after D; IO1 <= ‘0' after D; 
IOO <= ‘1' after D; --57 

when "01111111000101" => I05 <= '1' after D; I04 <= '1°* after D; 
IO03 <= ‘1' after D; I02 <= '0' after D; IO1 <= ‘1' after D; 
IO0 <= '0' after D; --58 

when “J1ITII1110O001G1* =s"° 105 <= “1° after DD; 104 <= *1' -atter D; 
I03 <= '1' after D; I02 <= '0' after D; IOl1 <= ‘1' after D; 
IOoO <= '1' after D; --59 

when °11111111000110" => I05 <= ‘1' after D; I04 <= ‘1' after D; 
1O3 <=" |“ -atter D> 102) <= '1* atter D;- I0l. <= "0° after D; 
IoO0 <= '0' after D; --60 

when °01111111000110" => IOS <= '1' after D; I04 <= ‘'1' after D; 
IO03 <= '1' after D; I02 <= '1' after D; IO1 <= '0' after D; 
Io0 <= '1' after D; --61 

when °01111111000111" => I0O5 <= ‘1' after D; I04 <= '1' after D; 
TOs <= "1watter D> 102 <= “L* atter D; I10l-<=:"'1° atter -D; 
I0o0 <= '0' after D; --62 

when *%00111111000111" => IOS <= ‘1' after D; I04 <= '1' after D; 
I03 <= ‘l' after D; IO2 <= '1' after D; IO1l <= '1' after D; 
Io0 <= '1' after D; --63 

end case; 
end process Pl; 
end EEPROM_2A; 


ee cm ee ee eee eee eee ee ew mm wm we oe ee ee i i ee 


--Test Bench Module for: EEPROM_2 
-- Note: comments beginning with WIZ should be left intact. 
-- Other comments are informational only. 


-- Multisim VHDL software version: 3.75a 


library ieee; 
use ieee.std_logic_1164.all; 


-- use ieee.numeric_std.all; -- Note: uncomment this if you use 
-- IEEE standard signed or unsigned types. 
-- use ieee.std_logic_arith.all; -- Note: uncomment/modify this if you use 


-- Synopsys signed or unsigned types. 
use std.textio.all; 
use work.all; 
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entity TESTBNCH is 
end TESTBNCH; 


architecture stimulus of TESTBNCH is 
component EEPROM_2 is 
Dore o( 

AO: in BIT; 
Ad? an BLT; 
A2: in BIT: 
AS? 1n- SIE: 
A4: in BIT; 
A5: in SIT; 
A6: in BIT; 
Ag+ an Bit; 
AS: in BIT; 
A9: in BIT; 
AlO: in BIT: 
All: 2n- BIT; 
Al2: in BIT; 
Al3: an BIT: 
LO0 = Out BIT; 
LOL: Out BIT; 
TO2: Out BIT: 
103% Out BIT; 
104: ouce BLT: 
1O5:* ouc BIT 


‘3 


end component; 
constant PERIOD: time := 50 ns; 
-- Top level signals go here... 
Signal AQ: BIT; 
Signal Al: BIT; 
Signal A2: BIT; 
Signal A3: BIT; 
Signal A4: BIT; 
Signal AS: BIT; 
Signal A6: BIT; 
Signal A7: BIT; 
Signal A&8: BIT; 
Signal A9: BIT; 
Signal A10: BIT; 
Signal All: BIT; 
Signal Al2: BIT; 
Signal Al13: BIT; 
Signal IOO0O: BIT; 
Signal IO1: BIT; 
Signal 102: BIT; 
signal IO3: BIT; 
Signal I04: BIT; 
signal IO5: BIT 


Signal done: boolean := false; 


begin 
DUT: EEPROM_2 port map ( 
AO, 
Al, 
A2, 
A3, 
A4, 
AS, 
A6, 
A7, 
Ag, 
AY, 
Alo, 


) 


STIMULUS]: process 


-- Sequential stimulus goes here... 


begin 

AQ <= 'Q', 
Al <=a"l1"> 
AZ <= "OC", 
As <=) °2" 3 
Aa <=) "Cs 
A5 <= '0'; 
AG <= *1', 
AD <6] i; 
A8 <= 'l', 
AS <= ‘l', 
ALQY<= "1". 
Ali tea 1s, 
Ril2gs<= Oo! 
Al3 <= '0O'; 
wait; 


Ald, 
Al2, 
ALS; 
Io0, 
101; 
LOZ; 
1G 
I04, 
I05 


t 


Sample stimulus... 


end stimulus; 


al 
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after 


after 


after 
after 
after 
after 


after 100 ns; 
after 100 ns; 


100 


100 


100 
100 
100 
100 


ns; 


ns; 


ns; 
ns; 
ns; 
ns; 


-- Suspend simulation 
end process STIMULUS1; 


oS 
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APPENDIX B. MATLAB SIMULATION AND DATA 
EVALUATION CODE 


$RSNS8 WITH SHIFTING 


% LT Nathan York, 15 December 2000 

% Generates the simulated folding waveform. Two element arrays 
assumed. 

% Can accept normalized mixer output for comparison to simulated 
waveform. 

% Plots folding waveforms for simulated and measured data. 

% m = modulus; d = spacing (m); f£ = frequency 

clear 

N=2; 

f=8.0e9; c=3e8; wavl=c/f; 

k=2*pi/wavl; 

rad=pi/180; 

sf=sqrt(3)/2 % sf is the scale factor 


% ffd is forced phase difference to align signal to required zero 
crossing 
ffd=45; toriginal is 45! 


% number of moduli in array design and their values 
mm=[8 17); nm=length(mm) ; 


disp('moduli are: '),disp(mm) 
% specify a modulus to plot 
m=S-: 


% Mis Dynamic Range 
M=64; 

% nf is number of folds 
nf=M/ (2*m*N) ; 


% determines required element spacing 
d=nf*wavl/sf/2; 


%$ start, stop, increment angles 
Starct=-90; stop=90- ine=0.1- 
N=floor((stop-start) /inc) +1; 


Gorm m=: he 
thd(n)=start+(n-1)*ince; thr=thd(n) *rad; 
arg(n) =Kk*a*sin (thr); 
%$ mixer output 
mx (n) =cos (arg(n)+(ffd)*(pi/180) ); 
end 
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load rsnsi4.dat 
boths=(|[thd*  mx*] > 
Save simu8.dat both8 -ascii 


figure (1) 

hold 

DLoetend ms, 5°): 
pPlectrsnsi4 (471) yrsnsi4 ys) 


hold off 

xlabel (‘Angle of Arrival'),ylabel('Normalized Mixer Output'); 
grid 

xlabel('Angle of Arrival (degrees) ') 

ylabel('Normalized Mixer Output') 

title(i*Modulus 8, with fid = “mumZstr(fid) ]) 
legend('Simulated', 'Measured',0Q) 

axis([{-90,90,-1,1]) 

orient landscape 


temp=rsnsi14(:,3); 
temp1l=rot90 (temp) ; 
left=templ-mx; 
zero=0; 


figure (2) 

hold 

Dise(Ends noc, -b" > 
pleottvsnsi4(:, 1) rsnsi4(. 3) 5 
Dlott{thd, vert, 'g'); 

DiSot( tnd, zero, 'k"): 

xlabel('Angle of Arrival (degrees) ') 
ylabel('Normalized Mixer Output') 
legend('Simulated', 'Measured', 'Difference',0) 
axis([{-90,90,-1,1]}) 

Grid 

title('Modulus 8, Measured values - Expected values') 


%RSNS17 WITH SCALING 


% LT Nathan York, 15 December 2000 

% Generates the simulated folding waveform. Two element arrays 
assumed. 

% Can accept normalized mixer output for comparison to simulated 
waveform. 

% Plots folding waveforms for simulated and meaSured data. 

%m = modulus; d = spacing (m); f = frequency 


clear 
N=2; 
f=8e9; c=3e8; wavl=c/f; 
sf=.8746; %sqrt(3)/2 
k=2*pi/wavl; 
rad=pi/i80; 
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% ffd is forced phase difference to align signal to required zero 
crossing 
ffd=90 %original is 90! 


m=17 ; 


% M is Dynamic Range 

M=64; 

% nf is number of folds 

n£=M/ (2*m*N) ; 

% calculation of required element spacing 
d=nf*wavl/sf/2; 


% start, stop, increment angles 
Sstart=-90; stop=90; inc=0.1; 
N=floor( (stop-start) /inc)+1; 


£Or- n= eN- 
thditn) =Stane+(m-1)*ince;> thr=thd(n) *rad; 
arg(n) =kK*d*sintehr); 
% mixer output 
mx (n)=cos(arg(n)+(f£f£d) * (pi/180)); 


end 


load rsnsll.dat; 
load rsnsl14.dat; 
load PAT9.dat; 

load PAT10.dat; 
load patl0a.dat 


Poem 7= (tha mx ae 
save simul7.dat bothl7 -ascii; 


figure (3) 
hold 
DIOELENG, mx, 'D )4 


% LABVIEW Output has mod 17 mixer output in second colum. 

% Normalization technique puts normalized mixer output in third 
column. 

elot(rsnsl4(:34),rsnsl4(:.6),'r') 


hold off 

xlabel('Angle of Arrival (degrees)'),ylabel('Normalized Mixer Output') ; 
grid 

title([{'Modulus 17, with ffd = ',num2str(ffd)]) 


legend('Simulated', 'Measured',0) 
xlabel('Angle of Arrival (degrees) ') 
ylabel('Normalized Mixer Output') 
axis ([-90,90,-1,1)) 

orient landscape 


oi 


temp=rsnsl14(:,6); 
templ=rot90(temp) ; 
left=templ-mx; 
zero=0; 


figure (4) 

hold 

Dilot (end mx, bb") 

plet (rsnsl4(:,4),rsnsl4(:,6),'‘r*} 
plOEl thd plerr, ¢')} 

plot (CEha- zero; 'k') 
legend('Simulated', 'Measured', 'Difference',0) 
Grid 

xlabel('Angle of Arrival (degrees) ') 
ylabel('Normalized Mixer Output') 
axis([(-90,90,-1,1])) 

title('Modulus 17, Measured values - Expected values') 


$BinPop.m 
%$ LT Nathan York 15 December 2000 


%$ step size is 0.1 degrees based on minimum stepper motor resolution 
% simulates the entire transfer function 

% 1- determines number of active comparators at each step for both 
modulus 

% 2- determine bin number 

% 3- determine angle of arrival 


% Also calculates bin population and quantization error. 


% Load either the simulated data (simul7.dat, simu8.dat) 

% or the experimental data - SEE BELOW (ALSO SEE LINE 155-156) 

% Uses predicted thresholds calculated from cosine of equally spaced 
angles 

Bek RK KK ek kkk kk kk kk kk kk kk ke kk 

format long; 

clear; 

sum_error_squared=0; Excess_RE=0; 

Sfi=sqrt(3) /2; 


$load simul7.dat 
$load simu8.dat 
$df1=simul7; 


load rsnsl14.dat 
afl=(rsnsl4(:,4) rsnsl4(:,6)]; 


ERKKKKKKKKKEKKKEKEKKKKK KEE KK KK KK KK KK 


point_num=1801; 
Sum_AE=0; 


$comparators mod 17 
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fOr Aa We porte num 


e2mi? (io) h)) 
Cum (ie ez it) 
em) 7 (ie 3 
© Fmd a) 
oven ey Glen Eo) |e) 
c_m17(i,[6])= 
eran le Cael 1) 
ecomlL7 (ie) 
elm? (25.091) ) =e; 
Coml? (i (10) —0 - 
© mi7 (2, [11])=0; 
Goma) (1, [1 21) =0; 
Com 7 (i (031) =0- 
eomL7 (25114) 7=0- 
mel (in [eo =0- 
Gi (ie (6) =0: 
emi? (a, [7 1-0; 
Gum. 7 (3a (ie | =o: 


end 


S%SIMULATION THRESHOLDS 
tl = cos(33*pi/34); 
t2. 3= cos(3i*pi/3s4)- 
ts. = Cost 29" p17 34) 
Lo =Cos (27 =p / 34); 
5 “= .«Cos (25° p17 34); 
t6 = coes(Z23*pi/34); 


t7 2 = Cos(Z1i*pi/34)- 
Go -="CcOs (19 * pic) 3 
tO = COS (17*p1/34)- 


SLOs=. cost ls *pi7 32); 
Cli, = Costissp1/s2) 
t12 = cos(11*pi/34); 
ti3: = Cos (9"*p1/ 34); 


t14 = cos(7*pi/34); 
G15) = .€oOs (97 p17 54):- 
tl6o- = e€os(3*p1/34); 
tl? ‘= se0s(p2/34); 


Boo i=) >polnt. num 
% all the data points 


% cl is the bottom comparator 
% tl is the threshold voltage for cl etc. 


gE GE (awl 2) eet 


Com GA eas 
else 
Com? Gel Wee= 00> 


end 
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Le ob e215) See 2 
Gsm? (ila) = ake 
else 
esl (as = 
end 


ph emesis (21 fe) Bm | 
eoml7 (i. (4) ) a2: 
else 
eoml7 (2,(4)) >= 70- 
end 


HE (oleae ah) ) > eS 
Caml7 (i (Sl) e= 1; 
else 
Sami (1715) ) <=40% 
end 


De Sta eae i. 2) aa Go 
Cami (i561) tae: 
else 
eomity (25 (6) ) =20> 
end 


ee (Cit al 2 eit S 
Smee) 7 lee 
else 
Camis )) «= .0- 
end 


re (ab Ca 2 ye te 
oom (2, [Si jv= 1 
else 
emi s41.1S)) = 0; 
end 


Lt dt lat) 3 
Gami7 i. 191) = 1 
else 
Cami 779) -=- 0; 
end 


ny eee 6 baad Ga a (ae 
e2m17 41 (| 20-1 
else 
Gorm 7 (22120 pa: 
end 


Dore e Lele) ) et 
erm 7 |) a 
else 
eomly (1 j-[124) 
end 
be ae a, (2) bl 


1 
© 
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Comey Ae elon ae 
else 

Com (1. bl 2)) c=. 0; 
end 


ier Ati 2) eZ 


Canis (ils) ee 
else 

Camis = Oe 
end 


ay ag (ob 2B se Bee oa he 


Cun (aie (a4 je7= is 
else 

Cum fees) = “0 
end 


PES ECUE 1 UA ez eat he 


Cums Wee ora a= cy 
else 

Gumi (a1) = -0- 
end 


PE eK Ge at) eS 


‘odmesgy il lot Glued) [ile o\0 (> ame ame Oe 
else 

Cop tas G)) = 20; 
end 


DP, Carr, 2) jet 6 
Caml tivity hy cers 
else 
Go mily kay) 
end 


i) 
© 


de (Od Cae) ie Ly 


Caml | 8 eee 
else 

Comlsy(1, (181) 2= 0: 
end 


end 
encoding the thermometer code into a decimal number. 
$rsns_m17 contains the number of comparators on at each time 
$step 


Lor j=l: point num 


Recoml 7 (34 02 ))):5 
B=c_m17(3,(3]); 
C=e ml7(3,([41]); 
D=c_m17(3,([5]); 
E=c_m17(3,(6]); 
F=c_m17(3,([7]); 
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G=eEr7ml 7 (7 el) s 
Becuml (5, Lol); 
Lee mi? 104) 
J=erml yy 
Keer ml ey tle ee 
PHeum jul sa es 
M=camd) (34-1441). 
N=Goml7 (44S 


O=caml 717, (161); 
P=Cami7 (j, 1271) 
Q=e.ml7(7,(18)])-; 


Z=A+B+C+D+E+F+G+H+1I+J+K+L+M+N+0+P+Q; 
rsns_m17 (j)=2; 

end 

xarray=[-90:180/1800:90]; 

figure (6) 

plot (xarray,rsns_m17(:)),grid 

xlabel('Angle of Arrival (degrees) ') 

ylabel('Number of Comparators') 


% title('Number of Comparators on each time step for mod 17') 
axis([{-90 90 -1 17]) 


BKK KKK KKK KKK KKK KK KKK KEK KK KKK KK KK 


clear dfl 
$dfl=simu8; 


dfil=(rsnsil4(:,1) rsnsl4(:,3) 1]; 


Bxexxekkekeeeae keke ke KK KEK KK KKK KK KEK KKK 


$comparators mod 8 


fOr 2i=1 *point- num 


tl = cos(15*pi/16); 


t2 = cos(13*pi/16); 
t3 = cos(11*pi/16); 
t4 = cos(9*pi/16); 
C5 = ¢os(7*pi7 16); 
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C6)=" €COS(5*p17 16); 
G7 = \€0S(3"*pi/16) > 
ESu= “COS(p1/ 16) 


for i=l1:point_num 


% all the data points 
% cl is the bottom comparator 
% tl is the threshold voltage for cl etc 


5s Seay To al Lt Gs Ww 0 2 > a ob 
come (2, [2)j2 = 1; 
else 
Gamevas [25> 
end 


tl 
) 


Pt edhe 2) ye 2 
Game eastouees = 
else 
come Ci, [3)) 
end 


If 
o 


TE (didi. 12))) > cs 
eome@wttl)) = 1: 
else 
e_m8(i,[4)) 
end 


tl 
>) 


aE (d£l(i,(2)))> t4 


emo clo) ee 
else 

Come (ie Ses Oo 
end 


Dah ce oh cael is oy (Ped (oo eas 
ems(2, (6])--= 1; 
else 
c_m8 (i, [6]) 
end 


ll 
>) 


Vf (delice) ))= 6 
Come (i ty |) =) LL 
else 
Camst 7 1) 
end 


il 
>) 


Le tat biti ee bk 7 
ene (oles Ls 
else 
@lmg ( 12,03) )* == 0: 
end 


if (dhl alZh) Se cs 
Cano o l= 2 
else 
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Gams (i (9 e— a0: 
end 


end 


encoding the thermometer code into a decimal number. 
$rsns_m8 contains the number of comparators on at each time 
$step 


fOr j=12 point. num 


A=c_m8(j,(2]); 
B=c_m8(j,[3]); 
C=c_m8(j,[4]); 
Dee. ms), bol); 
E=c_m8(j,[6]); 
F=c_m8(j, [7] i 
G=c_m8(j,[8]); 
H=c_m8(j,[9]); 
Z=A+B+C+D+E+F+G+H; 


rsns ms ())=2- 
end 


figure (7) 

plOe (xearray..sSns= mst)’ ©"). grid 

xlabel('Angle of Arrival (degrees) ') 

ylabel('Number of Comparators') 

% title('Number of Comparators on each time step for m8') 
axis([-90 90 -1 8])) 


figure (8) 

hold on 

plot (xarray,rsns_ml17(:),'b') 

plot (xarray,rsns_m8(:),'r') 

Noda hots 

xlabel('Angle of Arrival (degrees) ') 
ylabel('Number of Comparators') 
legend('Mod 8','Mod 17',0) 
$title('Number of Comparators on each time step for mod 17 & mod 8') 
axis({-90 90 0 17)) 

orient landscape 


% integers for trans_funct file 
bin=-1; bin_ctr=1; start_th=-90; end_th=0; sum_error_squared=0; 


% Bin mapping for thermometer code. doa is the bin number. 
For s=1:point num 


1f (rsns_ml17(s)==10 & rsns_m8(s) ==6) 
doa(s) = 0; 


elseif (rsns_m17(s)==10 & rsns_m8(s)==5) 
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doa(s) = 1; 

elseif (rsns_m17(s)<==11 
doa(s) = 2; 

elseif (rsns_ml17(s)==11 
doa(s) = 3; 

elseif (rsns_ml17(s)==12 
doa(s) =4; 

elseif (rsns_ml17(s)==12 
doa(s) = 5; 

elseif (rsns_m17(s)==13 
doa(s) = 6; 

elseif (rsns_ml17(s) ==13 
doa(s) = 7; 

elseif (rsns_m17(s)==14 
doa(s) = 8; 

elseif (rsns_ml7(s)==14 
doa(s) = 9; 

elseif (rsns_m17(s)==15 
doa(s) = 10; 

elseif (rsns_ml17(s)==15 
doa(s) = 11; 

elseif (rsns_ml17(s) ==16 
doa(s) = 12; 

elseif (rsns_ml17(s)==16 
doa(s) = 13; 

elseif (rsns_m17(s) ==17 
doa(s) = 14; 

elseif (rsns_ml17(s) ==17 
doa(s) = 15; 

elseif (rsns_ml17(s)==16 
doa(s) = 16; 

elseif (rsns_m17(s)==16 


doai(s) = 17; 


elseif (rsns_ml17(s)==15 
doa(s) = 18; 


elseif (rsns_ml17(s)==15 
doa(s) = 19; 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


& 


rsns_m8(s) ==5) 


rsns_m8 (s) ==4) 


rsns_m8(s)==4) 


rsns_m8(s) ==3) 


rsns_m8 (s) ==3) 


rsns_m8 (s) ==2) 


rsns_m8(s) ==2) 


rsns mo(sy==1) 


rsns_m8(s) ==1) 


rsns_m8 (s) ==0) 


rsns_m8(s) ==0) 


ESnSenc (Ss )e= 1} 


rsns_m8(s) ==1) 


rsns_m8 (s) ==2) 


rsns_m8 (s) ==2) 


rsns_m8 (s) ==3) 


rsns_m8(s)==3) 


rsns_m8 (s) ==4) 
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elseif (rsns_ml17(s)==14 
doa(s) = 20; 

elseif (rsns_ml17(s)==14 
doa(s) = 21; 

elseif (rsns_m17(s)==13 
Gdoa(s) = 22; 

elseif (rsns_ml17(s)==13 
doa(s) =). 23; 

elseif (rsns_ml17(s)==12 
doa(s) = 24; 

elseif (rsns_m17(s)==12 
deats) —=—25. 

elseif (rsns_ml17(s)==11 
doa(s) = 26; 

elseif (rsns_m17(s)==11 
doa(s) =" 27; 

elseif (rsns_ml17(s)==10 


doa(s) 28; 


& 


& 


elseif (rsns_m17(s)==10 & 


doa(s) 29; 
elseif (rsns_ml17(s)==9 
doa(s) 20: 


(rsns_ml17(s) == 
31> 


elseif 
doa(s) 


elseif (rsns_ml17(s)== 
doa(s) 32 


elseif (rsns_m17(s)== 
doa(s) 335 


elseif (rsns_ml17(s)== 
doa(s) 34; 


elseif (rsns_m17(s)== 
doa(s) 35: 


elseif (rsns_m17(s)== 
doa(s) 36; 


elseif (rsns_m17(s)== 
doa(s) Sa: 


& 


& 


& 


& 


& 


& 


& 


rsns_m8 (s) ==4) 


rsns_m8 (s)==5) 


rsns_m8(s)==5) 


rsns_m8 (s) ==6) 


rsns_m8 (s) ==6) 


rsns_m8 (Ss) ==7) 


rsns_m8 (s) ==7) 


rsns_m8 (s) ==8) 


rsns_m8 (s) ==8) 


7) 


rsns_m8(s)= 


rsns_m8 (s) ==7) 


rsns_m8 (s) ==6) 


rsns_m8 (s) ==6) 


rsns_m8 (s) ==5) 


rsns_m8 (s)==5) 


rsns_m8 (s)==4) 


rsns_m8 (s) ==4) 


rsns_m8 (s) ==3) 


106 


elseif 


(rsns_mil17 (s) == 


doa(s) = 38; 


elseif 


(rsns_m17(s)== 


doa(s) = 39; 


elseif 


(rsns_ml17(s) == 


doa(s) = 40; 


elseif 


(rsns_m17(s) == 


doa(s) = 41; 


elseif 


(rsns_m17(s)== 


doa(s) = 42; 


elseif 


(rsns_ml17(s)== 


doa(s) = 43; 


elseif 


(rsns_ml17(s) == 


doa(s) = 44; 


elseif 


(rsns_m17 (s) == 


doa(s) = 45; 


elseif 


(rsns_m17(s) == 


doa(s) = 46; 


elseif 


(rsns_m17(s) == 


doats) = 47; 


elseif (rsns_ml17(s)==0 


doa(s) = 48; 


elseif 


(rsns_ml17(s) == 


doa(s) = 49; 


elseif 


(rsns_ml17(s)== 


doa(s) = 50; 


elseif 


(rsns_m17(s)== 


doa(s) = 51; 


elseif 


(rsns_ml17 (s) == 


doa(s) = 52; 


elseif 


(rsns_ml17(s) == 


doa(s) = 53; 


elseif 


(rsns_mi17(s) == 


doa(s) = 54; 


elseif 


(rsns_m17(s) == 


doa(s) = 55; 


elseif 


(rsns_ml17(s)== 


rsns_m8(s) ==3) 


rsns_m8(s) ==2) 


rsns_m8 (s) ==2) 


rsns_m8(s) ==1) 


rsns_m8(s)==1) 


rsns_m8(s) ==0) 


rsns_m8 (s) ==0) 


rsns_m8 (s) ==1) 


rsns_m8 (s) ==1) 


rsns_m8(s) ==2) 


rsns_m8 (s) ==2) 


rsns_m8 (s) ==3) 


rsns_m8 (s) ==3) 


rsns_m8 (s) ==4) 


rsns_m8(s) ==4) 


rsns_m8 (s) ==5) 


rsns_m8(s)==5) 


rsns_m8 (s) ==6) 


rsns_m8 (s) ==6) 
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doats) =.56; 


elseif (rsns_mi7(s)==4 & rsns_m8(s)==7) 


Goats) h=<57: 

elseif (rsns_mi7(s)==5 & rsns_m8(s) ==7) 
doa(s) = 58; 

elseif (rsns_ml17(s)==5 & rsns_m8(s) ==8) 
doa(s) = 59; 

elseif (rsns_mi7(s)==6 & rsns_m8(s)==8) 
doa(s) = 60; 

elseif (rsns_mi7(s)==6 & rsns_m8 (Ss) ==7) 
doa(s) = 61; 


elseif (rsns_ml17(s)==7 & rsns_m8 (Ss) ==7) 


Goats)” = 62; 
elseif (rsns_mi7(s)==7 & rsns_m8(s) ==6) 
doa(s) = 63; 
else 
doa(s) = 100; % if the vector does not map, 100 is the error 
code 
end 


% Uses the formula asin{[(#-34)/64] to find bin center 
gamma=doa(s)+1; 
bin_center (s)=2*asin((gamma-32) *sf/64) *180/3.1415; 


Bole 
1f doa(s)==100 
Din .center (s)=0-: % if the vector did not map, reports as 
broadside 
end 


angle_error(s) =bin_center(s)-(s-900)/10; 


Calculation of RMS Reporting Error 
$if statement eliminates excessive errors and subs with mean error. 
if abs(angle_error(s)) <= 15 
Sum_AE = Sum_AE + angle_error(s) *angle_error(s); 
Excess_RE=Excess_RE+1; 


end 
end 
Z2=0. 63: 
for -l=1 264 
bin count (1) =0- 
end 
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ave penile: Aa ilaallesy Ola 
if doa(lb)==100 
De= 1; 
else 
bc=doa (lb) +1; 
end 
bine Count (be)}=bin- count (bc) +1: 
end 


figure (9) 
hold 
Plec(2, bin count, 'b'), axis (10,63:0;60)]) ,<grid 
NG laAOrtE 
xlabel('Bin Number' ) 
ylabel('Bin Population') 


figure (10) 
plot (xarray,doa,'b'),axis([-90,90,0,65}) 
xlabel('Angle of Arrival (degrees) ') 
ylabel('Bin Number') 
grid 


figure(11) 
plot (xarray,bin_center,'b'),grid 
xlabel('Angle of Arrival (degrees) ') 
ylabel (‘Reported Angle of Arrival (degrees) ') 


figure (12) 
plot (xarray,angle_error, 'b'),grid 
axis([-60,60,-10,10]); 
ylabel('Reporting Error (degrees) ') 
xlabel('Angle of Arrival (degrees) ') 


% end of program. 


% phase_error.m 


% Determines the phase error between a predicted and measured waveform 
% 2 channel input 
% Written by Dave Wickersham 


clear all; 


load rsnsl14.dat 
theta=([rsns14(:,1)]j; 


eve—([rsnsl4(-,3) 1; 
avilJ=(rsnsia(:. 6) |. 
md=length(dv8) ; 


% constants and design variables 
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N-2e 

f=8e9; c=3e8; wavl=c/f; 

St=2/saqrt (3); 

k=2*pi/wavl; 

M=64; 

nf=M/ (2*N);% nf is number of folds (not really, no modulus in this 
equation) 

dw=nf/2; % new design for mixer ouputs 


FirstMin=0;FirstMax=0; SecMin=0; SecMax=0 ; Max17=0;Min17=0; 
rad=pi/180; 
thet_r=theta*rad; 


ffd8=45; $45 
ff£d17=90; $90 


e8_sq=0; e17_sq=0; cnt8=0; cnti7=0; 
da8=dw*wavl*sf/8: 
di7=dw*wavl*sf/17; 


subroutine to find local max and min 
% works because I know about where the max and min will fall 
% used to get the right quadrant for the arccos 
%Mod 8 
for i=(md-1)/6:(md-1)/2 
if dv8(i)<FirstMin 
FirstMin=dvs (1); 
Index1Min=i; 
end 
end 


£O6r W=(ma—-1) 7/3: (ma-1)/ 2 
if dv8(1)>FirstMax 
FirstMax=dv8 (1); 
Index1Max=1; 
end 
end 


for i=(md=1)/2:2* (md-1) 73 
if dv8(i)<SecMin 
SecMin=dv8 (i); 
Index2Min=1; 
end 
end 


FOr d=2*(ma—-1)/3-3* (ma-1)79 
if dv8(i)>SecMax 
SecMax=dv8 (1) ; 
Index2Max=1; 
end 
end 


Mod 17 
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fOr 1=(ma-1)7225* (mad-1)/6 
re nevi?s (ijeMini? 
Mini? =dv1l7 (1)? 
Index17Min=1; 
end 
end 


for 1=(md-1)/6: (md-1)/2 
if dv17(1)>Max17 
Max17=dv17 (1); 
Index17Max=1; 
end 
end 


tor n=1- ma 

delt8 (n) =k*d8*sin(thet_r(n))+(f£d8) *rad; $delt8 is the simulated 
total angle in rad 

unwrap (delt8(n) ) ; 

while delt8(n)>pi 

% delt8 (n) =delt8(n)-2*pi; 

$end 

while delt8(n)<-pi1 

% delt8(n) =delt8(n)+2*pi; 

$end 


am8 (n) =acos (dv8(n) ); $dm8 1s measured data's total angle in rad 
% unwrap (dm8(n) ; 
1f (n>=IndexlMin) & (n<=Index1Max) 
dm8 (n) =-dm8 (n) ; 
end 
if (n>=Index2Min) & (n<=Index2Max) 
adm8 (n) =-dm8 (n) ; 
end 


phi8(n)=(delt8(n)-dm8(n)); 
% wumnwrap phi8(n); 

while phi8(n)>pi 

phi8 (n)=phi8(n)-2*pi; 
end 
while phi8(n)<-pi 

phi8 (n) =phi8(n)+2*pi; 
end 


delt17 (n) =k*d17*sin(thet_ri(n))+(f££d17) *rad; $delti7 is the 
Simulated total angle in rad 
% unwrap (delti17(n)); 
while delt17(n)>pi 
delt17 (n)=delt17(n)-2*pi; 
end 
while delt17(n)<-pi 
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deéltl7 (na) =dele17 (n) +2*01- 
end 


dmi7 (n) =acos (dv17(n)); ¢dm17 is measured data's total angle in 
rad 


if n>=Index17Min 
admi7(n)=-dm17 (n); 
end 


if n<=Index17Max 
adm1i7(n)=-dm17(n); 
end 


phil? (n) =(delti7(n)-dm17(n)); 
while phil7(n)>pi 

Phil? (n)=phil7 (n) =2* pi; 
end 
while phil7(n)<-pi 

Phil? (nh) =p 7 in) +2 "pi; 
end 


end 


phi8=phi8/rad; 
Dh if =p vad; 


Eon 1 => iad $phi contributes to RMS error 
only if it is real 

e8_sgq = e8_sq + (real (phi8(i))*real(phi8(i))); 

el7_sq = el7_sq+ (real(phil7(i))*real(phil7(i))); 
end 


rms8=sqrt (e8_sq) /(md-1) ; calculate RMS phase error for each 
channel 
rmsl17=sgqrt(e1l7_sq) /(md-1) ; 


disp(['RMS Phase Error of Mod 8 is ',num2str(rms8),' ']) 
disp({'RMS Phase Error of Mod 17 is ',num2str(rmsi7),' ')) 
figure(5) 

hold 

SDlOC( D126: 4) ,D12(2,7 5), —-t  ) 


plot (theta,phi8,'b'); 

plot (theta,phil7,'r'); 

axis ([-90,90,-90,90])) 

Holavorcer 

legend({'Mod 8', 'Modi7',0) 

xlabel('Angle of Arrival (degrees)'),ylabel('Phase Error (degrees) '); 
Grid 
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$title([{'Phase Error between Simulated and Measured Data, Pattern 32']) 
orient landscape 


ie 


THIS PAGE INTENTIONALLY LEFT BLANK 
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APPENDIX C. PRINTED CIRCUIT BOARD CONSTRUCTION 


The following is a list of steps that were used to create the printed circuit boards 
for this thesis. These steps can streamline the construction process for someone who is 
not familiar with the process. The instructions are only valid when using the Electronic 
Workbench Multisim software package and the PCB milling machine operated by the 
Physics Department at NPS. 


1. Start Multisim. Create and test the circuit you would like to build. Be sure to include 
all power sources and grounds. 

2. Select the Transfer to Ultiboard option in the Transfer pull down menu. This will 
automatically start the Ultiboard program and import a parts list and netlist that is 
used by Ultiboard. 

3. Setup the Ultiboard workspace for your design. To do this you will have to select a 
design class, board size, number of layers, and trace size. Use the defaults for all the 
other options. 

4. You will notice a the outlines for all the components at the top of the screen. You 
will have to manually place each of the components within the board outline. Each 
component will have a line ending in a circle emanating from the part's center. This 
is a placement, aid which will help you determine the optimum placement for each 
component. The shorter the lead the better. As you place components in the board 
outline you will also notice yellow lines connecting the pins of the selected 
component and others. These are your pin connections. This can help you determine 
the appropriate orientation of the components. 

5. Ensure that all the power connections are on the board. These are the only parts that 
are not included by the Multisim program. You will need to manually add the 
connectors to the board and manually add each of the pins into the correct netlist. 

6. Once all the parts are in their postions save the design file and transfer the layout to 


the Ultiroute program. This program will automatically route all the traces for your 


hy hes 


10. 


jae 


board. In most cases, the first attempt may no complete all the traces. If this 
happens, close the program and go back to the Ultiboard program. You will need to 
reconfigure the location of the parts for a second run. If this happens, ensure that you 
check your grid size. Ultiroute sets the gridsize to 12.5 mils by default even if the 
design file uses a different grid setting. Make sure that you change the grid setting 
back to its onginal, otherwise the design will not transfer to Ultiboard for a second 


run. You need to continue this process until the board is 100% routed. 


. After the board is routed, transfer it back to the Ultiboard program. Use this program 


to check for continuity, power grids, pin connections and ground lines. One tip to 
make the construction easier is making all pin connections on the bottom layer. This 
will prevent you from having to solder underneath a component. 

Once this is complete, it is time to generate a Gerber file. This is a standard computer 
aided manufacturing (CAM) file type for all PCB milling machines. Under the File 
menu Select Post Processing. Here you will be prompted to the output device or file 
type. Select the Gerber RS-274-D (2.3) file type in the Photoplotter menu. Next 
double-click on each of the four layers to bring up the Plot Setting Parameters. 
Unselect the Extended Borders and Board Outline option. These option will cause 
the software used to run the milling machine to crash. The click on the Go! button. 
The Gerber files *.GO and *.G1 have now been created. The click on the Drill 
Targets button and add that information to the processed file. 

Copy the .GO, .G1 and .REP files to a 3 and a half inch disk and bring the disk down 
to the milling machine. 

The milling machine is running off a 486 computer using DOS 6.0. First read the 
-REP file and copy down all the apertures that are in the design file. These will have 
to be manually added to the Gerber file. 

Rename the .GO and .G1 files to TOP.GBR and BOT.GBR respectively. 


12. Copy these two files to the C:\CAM directory. 


13. 


Start the Gerber editor by typing CAM.EXE. 


14. Start a new design file and load the two .GRB files in it. 
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The next step is to edit the aperture list. The program will prompt you to do this. Just 
add the apertures from the .REP file to the list. 

Now it is time to edit the .GBR files themselves. Evaluate the design that was 
transferred from the Ultiboard program and make sure that all the traces and pads are 
present. At this point you will be able to modify the trace and pad sizes. I suggest 
that you use 20 mil traces and 80 by 80 mil square pads. These are the easiest to 
solder to. 

Once you are satisfied with the .GBR files, save your work and exit the program. 
Copy TOP.GBR and BOT.GBR to C:\PBOARD\INPUT. 

Rename TOP.GBR and BOT.GBR to L1.GBR and L2.GBR respectively. 

Using the EDIT program in DOS, open the L1 and L2 Gerber files and delete the G01 
from the third line and save. 

Type P1.BAT to start the outline converter. This program will convert the Gerber 
files to machine language that will control the milling machine. The outline converter 
may take up to two hours to convert a large design. 

Once the outline converter is done, type P3.BAT. This copies the output files to the 
CAM directory. 

Run the CAM program again. This time add LIOUT, L2OUT and DIOUT files to 
the original design file. You wll be able to check the new Gerber files against the 
originals to make sure the outline converter did its job properly. Most likely, the 
outline converter will miss a trace or pad. You will need to manually complete the 
trace in this program. If you need to change these files, make sure you go back and 
edit them to remove the G01 from the third line. Then copy the new files back to 
C:\PBOARD\OUTPUT. 

You are now ready to start milling. Type P2.BAT to start the milling program. The 
program will have prompts at the bottom of the screen to help you along. 

Place board material in milling machine. 

Start by setting the OPS Mode to Dnill. 

Select the layer DIOUT. 
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. Go to the Manual option to change out the drill bit and set drill depth. Use a 32 mil 


drill bit, and set the dnll depth so that the bit completely goes through the board 


material. Once this is set, quit and return to the main menu. 


29. Turn on the dnll and vacuum, and close the machine cover. The select Start. Follow 


30. 
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36. 


the prompts to start drilling all the holes. 

Once the holes are dnilled, it is time to cut the traces. To do this set the OPS Mode to 
mill and select the LIOUT layer. 

Go to manual to change the dmill bit to a milling bit. Set the depth of cut so that the 
channel is 10 mils wide. A couple of test cuts in the manual mode will be required to 


do this. 


. Go to Start and follow the prompts to start milling the traces. 
bi), 


Once the top layer is complete you need to cut the bottom layer. Follow the same 
procedure. The L2OUT layer needs to be selected as well as Mirror is ON selected. 
The bit and depth should be ok so go straight to Start and begin milling the bottom 
layer. 

Now that the board has been cut you need to check the board for milling defects. 
Start by sanding both sides of the board with a very fine grit sandpaper. 

Check for continuity and grounds. Do this using a multimeter. If a ground is 
detected, clean the channels with an exacto-knife. There is probably a small filing 
that was not completely removed by the drill bit that is causing the ground. 

When the board has been completely checked, it is time to stuff and solder it. Start 
by making the vias (through board connections). If there are not many vias, small 
pieces of wire can be soldered to both sides of the board. Another method is using 
small aluminum rivets. Put the rivets in the dll holes and turn the board upside 
down on a completely flat working surface. Use a punch to secure the rivet to the 
bottom of the board. Make sure that there is a good connection between the rivet and 
both sides of the board; otherwise you will have continuity problems. If you have 
trouble getting good contact with the rivets, a drop of solder on each side will correct 


the problem. 
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37. Now that the vias are complete, finish stuffing the board and soldering the pins and 
components into place. One hint: solder a few pins at a time and then check for 
grounds. It is very easy to accidentally jump the channel with the solder causing a 
ground. If check after a few pins it will be easier to find the mistake and correct it. 


38. The board is now ready to be used. 
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